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I. INTRODUCTION 


Q,3 


The LHC Run-I has confirmed the existence of a Higgs boson (h) [l|, whose properties 
are in agreement with those of the standard model (SM) within the uncertainties of the current 

n1 

data [3H8|- Thanks to the discovery of the Higgs boson, the SM was established as an effective 
theory to describe physics at the scale of electroweak symmetry breaking. In spite of the success 
of the SM, there are many motivations to consider new physics beyond the SM such as to solve 
the gauge hierarchy problem and to explain phenomena like neutrino oscillation, dark matter and 
baryon asymmetry of the Universe. There have been various new physics models proposed, some 
of which predict new particles at the electroweak to TeV scales. However, currently none of such 
new particles has been discovered yet. Their discovery is one of the main tasks of the LHC Run-H, 
which will start its operation in 2015. 

Even though the Higgs boson shows SM like properties, the Higgs sector can be extended 
from the minimal form with only an isospin doublet field. Indeed, there is no theoretical reason 
for the hypothesis of the minimal structure for the Higgs sector. Thus there are possibilities for 
extended Higgs sectors such as those with additional iso-singlets, doublets, and/or triplets. These 
extended Higgs sectors can also be consistent with all the current LHC data in some portions of 
their parameter space. 

Extended Higgs sectors are often introduced in various new physics models. Eor example, 
the Minimal Supersymmetric SM (MSSM) requires the Higgs sector with two doublet fields 9l.lid]. 
Multi Higgs structures are also studied in the context of additional CP violating phases [ll j and also 


realization of the strong first order phase transition 


I2I ]. both of which are required for successful 


electroweak baryogenesis [1^. Models with the Type-11 seesaw scenario are motivated to generate 


tiny neutrino masses by introducing a triplet field 


14l | . An additional singlet is required in the Higgs 


sector of the models with spontaneous breakdown of the U{1)b-l symmetry 15l4l7l|. which may be 
related to the mechanism of neutrino mass generation 181. Introduction of an additional unbroken 

nn 

symmetry into an extended Higgs sector, such as a discrete Z 2 symmetry [19|, [20(1 or a global U{1) 


symmetry 


2l|, can provide candidates of dark matter. Under the Z 2 or the global U{1) symmetry, 


if some of the scalar fields are assigned to be odd or to be charged, respectively, they cannot decay 
into a pair of SM particles so that the lightest one is stable. Such an unbroken symmetry can 
also be embedded into models with a radiative generation of neutrino masses [l^, 1221427] . where 


the existence of tiny neutrino masses and dark matter can be explained by the same origin of the 
symmetry. Therefore, a characteristic Higgs sector appears in each new physics model. 
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There are several important properties which characterize the structure of the Higgs sector. 
First of all, it is important to know the number of scalar multiplets and their representations. 
Second, does it respect new symmetries (global or discrete/exact or softly-broken)? Third, the 
mass of the second Higgs boson generally contains information of the new scale which does not 
appear in the SM. Fourth, the strength of the coupling constants among extra Higgs bosons provides 
information of the dynamics of the Higgs potential which is essentially important to understand 
nature of electroweak symmetry breaking. Finally, the decoupling property 2^ of extra Higgs 


bosons is closely connected to physics beyond the SM. Therefore, by future measurements of these 
properties, the Higgs sector can be reconstructed, and the direction of new physics beyond the SM 
can be determined. 

The direct search of extra Higgs bosons can provide a clear evidence to a non-minimal Higgs 


m 


30| 


sector. The current data accumulated from previous collider experiments such as LEP 
and Tevatron 3ll436t| have already given lower bounds for masses of the extra Higgs bosons. At 


the LHC Run-I, in spite of the discovery of a Higgs boson with the mass of 125 GeV, no extra 
Higgs boson has been found, and the parameter space for additional light Higgs bosons has been 


constrainec 
bosons 37 


to the considerable extent in regions with relatively smaller masses of the extra Higgs 


49l |. At the LHC Run-H, with the energy of 13-14 TeV and the integrated luminosity 


of 300 fb“^, wider regions of masses of the extra Higgs bosons will be surveyed. 

In addition to direct searches, new physics models beyond the SM have also been indirectly 
investigated by utilizing precision measurements of various physics observables such as the oblique 


parameters at LEP/SLC experiments 


0 - 


Elavour experiments have also been used to constrain 
the mass of charged Higgs bosons which appears in extended Higgs sectors M|, 1^. Now that the 
measured couplings of the Higgs boson h with the SM particles are consistent with the predic¬ 
tions in the SM within the uncertainties, it is time to consider fingerprinting of extended Higgs 
sectors 0, 54] by calculating radiative corrections to the predictions of those observables which 
will be measured with more precision at future experiments such as the LHC Run-H, the high lu- 
55l-l57t| with the integrated luminosity of 3000 fb“^ and future lepton colliders 


minosity (HL)-LHC 


like the International Linear Collider (ILC) l59(|. In new physics models with extended Higgs 


sectors, the coupling constants of h with the SM particles are generally predicted with deviations 
from the SM predictions due to field mixing and loop contributions of non-SM particles. Although 
no deviation has been found up to now in the Higgs boson couplings within the uncertainty of the 
current data, a deviation could be found in future experiments where more precise measurements 
will be attained. We then are able to indirectly obtain information of the second Higgs boson from 
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these deviations. Furthermore, a pattern of these deviations strongly depends on the structure of 
the Higgs sector, so that by comparing theoretical predictions of the Higgs couplings in various 
new physics models with future experimental data the shape of the Higgs sector can be determined 
indirectly. In order to compare the theory predictions to future precision data at the HL-LHC and 
also the ILC, where coupling constants are expected to be measured typically by a few percent 
or better accuracy, evaluations of the Higgs boson couplings including radiative corrections are 
inevitable. 

There are many studies for radiative corrections in extended Higgs sectors in the literature. 


Radiative corrections to the electroweak gauge boson two point functions (oblique correction^ 


tiave 


been studied in extended Higgs sectors in Refs. [6CH63l| . Loop induced vertices hgg [ 6 ^, /177 |65l-l7C| 
and / 1 Z 7 j^, 69-72| have been evaluated in extended Higgs sectors. Those to the Higgs boson 


couplings have been investigated in the two Higgs doublet model (THDM) in Refs. 73l-l76l| and in 
the Higgs triplet model in Refs. 0 ,S. 

In this paper, we study electroweak radiative corrections to the coupling constants of the 125 
GeV Higgs boson h in the THDM H with the softly-broken Z 2 symmetry 


symmetry, four types of Yukawa interactions 


8 fll |. Under the Z 2 


8ll-l84l| are possible depending on the assignment 


of the Z 2 charges into quarks and leptons. We investigate radiative corrections to the full set 
of Higgs boson couplings {hWW, hZZ, htt, hbb, hrr, hhh, / 177 , hZ'f and hgg) at the one-loop 
level in all types of the THDMs. We employ an improved on-shell renormalization scheme in our 
renormalization calculation where the gauge dependence in the calculation of the mixing angle in 
the previous studies is eliminated^. We then evaluate deviations in these coupling constants from 
the SM predictions under the constraint of current experimental data and theoretical bounds such 
as vacuum stability and perturbative unitarity. 

Furthermore, we investigate how we can extract information of the inner parameters such as the 
mass of the second Higgs boson and mixing angles when the scale factors are experimentally 
determined with the expected uncertainties at the HL-LHC and the ILC, where are the ratios 
of the measured couplings hXX from the SM predictions. Evaluating at the one-loop level in 
the THDMs, we discuss the possibility to measure properties of the Higgs sector using the future 
precision data by fingerprinting, and finally we determine the structure of the Higgs sector. 

This paper is organized as follows. In Sec. m we dehne the Lagrangian of THDMs, and give 
formulae for the Higgs boson masses and the Higgs boson couplings at the tree level. After that, 


^ According to Ref. [sst] . the gauge dependence exists in a renormalization of a mixing angle. 
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Z 2 charge 

Mixing factor 


‘I’l ‘1’2 Ql Ll Ur dR Br 


Type-I 

+ - + + --- 

cot (3 cot /3 cot j3 

Type-II 

-b - -b -b - -b -b 

cot P — tan P — tan P 

Type-X 

-b - -b -b - - -b 

cot P cot P — tan P 

Type-Y 

-b - -b -b - -b - 

cot/3 — tan/3 cot/3 


TABLE I: Charge assignment of the softly-broken Z 2 symmetry and the mixing factors in Yukawa interac¬ 
tions given in Eq. dH). 


we discuss constraints from vacuum stability and perturbative unitarity as the theoretical bounds. 
We then discuss the bounds from the electroweak oblique parameters, flavour experiments, direct 
searches of extra Higgs bosons at the LHC and the measurements of Higgs boson couplings at the 
LHC Run-I. In addition, we shortly summarize future prospects for extra Higgs boson searches 
and precision measurements of the Higgs boson h at the LHC Run-H, the HL-LHC and the ILC. 
In Sec. imi we explain renormalization in the electroweak sector, the Yukawa sector, and the Higgs 
sector in the THDMs. We also discuss the modified renormalization scheme. In Sec. m we give 
formulae of renormalized Higgs couplings and loop induced decay rates. We numerically estimate 
decoupling properties and non-decoupling effects of our one-loop calculations in the section. In 
Sec. |Vl we demonstrate how we can extract inner parameters by using future precision data. 
Discussions and conclusions are given in Sec. VI. 

II. TWO HIGGS DOUBLET MODELS 

A. Lagrangian 

In this section, we define the Lagrangian in the THDM with the softly-broken Z 2 symmetry, 
where the Higgs sector is composed of two isospin doublet scalar helds <I>i and ^ 2 - The charge 
assignment for the Z 2 symmetry is shown in Table U The following Lagrangian is modified from 
the SM: 

TtHDM = ^kin + Ty — V, (1) 

where Tkinj Ty and V are respectively the kinetic Lagrangian, the Yukawa Lagrangian and the 
scalar potential. Throughout the paper, we assume the CP invariance in the Higgs sector. 
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First, the kinetic Lagrangian is given by 


^kin = 


( 2 ) 


where is the covariant derivative: 

D, = d,- ^-gr^W^ - (3) 

with (a =1-3) and being the SU{2)i and U{1)y gauge bosons, respectively. The two 
doublet fields can be parameterized as 


^i = 


w- 


T /li T 


(* = 1 , 2 ), 


(4) 


where ui and V 2 are the vacuum expectation values (VEVs) for and $ 2 , which satisfy v = 
■\/v\^v\ = {\f2Gp)~^l‘^. The ratio of the two VEVs is defined as tan/3 = v^jvx. The mass 
eigenstates for the scalar bosons are obtained by the following orthogonal transformations as 








with i?(0) 


( cos Q — sin Q 
sin Q cos Q 



(5) 


where and G^ are the Nambu-Goldstone bosons absorbed by the longitudinal component of 
and Z, respectively. The mixing angle a is expressed in terms of the mass matrix elements 
for the CP-even scalar states as shown in Eqs. dUD-dii]). As the physical degrees of freedom, we 
have a pair of singly-charged Higgs boson , a CP-odd Higgs boson A and two CP-even Higgs 
bosons h and H. We define h as the observed Higgs boson with the mass of about 125 CeV. 

In terms of the mass eigenbasis of the Higgs fields, the interaction terms among the Higgs bosons 
and the weak gauge bosons are given by 

2 2 

Tkin =[sin(/3 - a)h + cos(/3 - a)i/] 

+ + 9<t,Yi,2ViV2 ( 6 ) 


where coefficients of the Scalar-Scalar-Cauge vertex and those of the Scalar-Scalar-Cauge- 

Cauge vertex g^^^^y^y^ are listed in Appendix A. 

Next, we discuss the Yukawa Lagrangian. The most general form under the Z 2 symmetry is 
given by 


—By —YuQl'^0'2^u'^ji + YdQ L^ddR + YeLh^eSR + h.c., (7) 
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where ^u,d,e are either ^ or $ 2 - Depending on the Z 2 charge assignment, there are four types 


of Yukawa interactions 


81 


82l |. which we call as Type-I, Type-II, Type-X and Type-Y 8^. The 


interaction terms are expressed in terms of the mass eigenstates of the Higgs bosons as 

-Y"‘= E + 


f=u^d,e 


^ - [VudU {mdid Pr - muiuPL) dH^ + me^e^PReH^ + h.c.] 


f f 

where K and are dehned by 


= sin(/3 — a) + cos(/3 — a), 
Ch = cos(/3 -a)-^f sin(/3 - a), 


( 8 ) 

(9) 

( 10 ) 


and C/ ia each type of Yukawa interactions are given in Table HI In Eq. (I8|) , represents the third 
component of the isospin of a fermion /; i.e., If = +1/2 (—1/2) for / = u (d, e). 

The Higgs potential under the softly-broken Z 2 symmetry and the CP invariance is given by 


V = + 77121 ^^ 2 !^ — m \{^\^2 + h.c.) 

+ + ^A 2 |$ 2 |^ + A3|$ip|^>2p + A4|^«l^>2p + ^As [($ 1 ^ 2 )' + h.c. 

The tadpole terms for hi and /12 are respectively calculated as 

Ti 

- = —m? + sin^ /3-(Ai cos^ /3 + Asin^ /3), 

V cos p 2 

T - 

-^ = —mo + cos^ /3-(A 2 sin^ (5 + Acos^ /?), 

u sin p 2 

where A = A 3 + A 4 + A 5 , and M describes the soft breaking scale of the Z 2 symmetry: 


( 11 ) 

( 12 ) 

(13) 

(14) 

sin /3 cos /3 

We note that can be taken to be both positive and negative values. By requiring the tree level 
tadpole conditions; i.e., Ti = T 2 = 0, mf and m^ can be eliminated in the Higgs potential. 

The squared masses of and A are calculated as 


= 


rrin 


m^± = - y(A4 + A5), = 

Those for the CP-even Higgs bosons and the mixing angle a are given by 

= cos^(a — fi)Mli + sin^(a — /3)M|2 + sin2(a — l3)Mi2, 
m\ = sin^(Q: — /3)Mfi + cos^(q: — /3 )M|2 — sin2(a — l3)Mi2, 

2MI2 


(15) 


(16) 

( 17 ) 












where Mf^ {i,j = 1,2) are the mass matrix elements for the CP-even scalar states in the basis of 
{hi,h2)R{/3): 

y2 _ 

Mil = cos^ P + X 2 sin'^ f3) -|- — Asin^ 2/3, (19) 

M 22 = M'^ + sin^ /3 cos^ /3(Ai -|- A 2 — 2A), (20) 

2 2 

Mf 2 = ^ sin 2;0(—Ai cos^ /3 + X 2 sin^ /?) + ^ sin 2/3 cos 2/3A. (21) 

Thus, ten parameters in the potential (^ 1^25 ^ 1.3 and Ai.s) can be described by the eight physical 
parameters rrih, a, /3, v and M^, and two tadpoles Ti and T 2 which are taken to 

be zero at the tree level. The quartic couplings A 1 -A 5 in the potential are then rewritten in terms 
of the physical parameters as 

Aiu^ = {m\ tan^ /3 -|- m\) sin^(/3 — a) -|- {m% -|- tan^ /3) cos ^(;8 — a) 

+ 2(m|^ — m|) sin(/3 — a) cos(/3 — a) tan /3 — tan^ /3, 

A 2 U^ = {rnjj cot^ (3 -|- m\) sin^(/3 — a) -|- {rnjj -|- ml cot^ /3) cos^(/3 — a) 

— 2{m\ — ml) sin(/3 — a) cos(/3 — a) tan f3 — cot^ /3, 

Asu^ = {m'jj — m|)[cos^(/3 — a) — sin^(/3 — a) -|- (tan (3 — cot /3) sin(/3 — a) cos(/3 — a)] 

-|- 2m‘^± — 

A4U^ = -|- m'ji — 2m^±, 

Asu^ = — m\. (22) 


We here dehne the so-called scaling factors to describe deviations in the Higgs boson couplings 


from the SM prediction as follows: 


Ky — 


„THDM 

i/hW 

slvv '' 


for H = Z, W, 


Kf = 


.,,THDM 

yhff 



\THDM 
_ '^hhh 
h — ■ 

'^hhh 


(23) 


where y^Yf and are the hVV, hff and hhh coupling constants in the SM, respectively, 


and those with THDM in the superscript are corresponding predictions in the THDM. The scaling 


factors for loop induced couplings can also be dehned by 


2 ^ r(h —)• 77 ) THDM 2 ^ r(h —Z 7 )thDM 2 ^ ^(h —)■ gg)THDM 

“ T{h 77)sm ’ ~ r(h ^ ^7)sM ’ ~ ^{h gg)sM 


where r(/i —)• Afy)sM and T{h —>• Xy)THDM are respectively the decay rates of the h —)• XY mode 








in the SM and in the THDM. At the tree level, the scaling factors are given by 


Ky = sin(/3 — a), 

i^f = ~ sin(/3 — a) + cos(/3 — a), 

K/j = sin(/3 — a) - ^^—— sin(/3 — a) cos^(/3 — a) 

M'^ — m? n , 

™ cos'^(/3 — a)(cot /3 — tan/3). 


ml 


(25) 

(26) 

(27) 


We can see that all the scaling 
this limit as the SM-like limit 


factors become unity when sin(/3 — a) = 1 is taken, so that we call 


86 |. 


It is convenient to introduce a parameter x defined as 


TT 


- - a), 


(28) 


where a; —>• 0 corresponds to the SM-like limit. We note that in the MSSM, the sign of x is 
determined to be negative due to supersymmetric relations [l^. Because the current LHC data 
suggest that the observed Higgs boson is SM-like, the case with |x| ^ 1 describes such a situation. 
In this case, we obtain 

^2 


Ky — 1 — + O(x^), 

x^ 

Kf = 1 + X — + 0{x^), 

/3 2M2 

1 “ 


ml 


x^ + 0{x-^). 


(29) 

(30) 

(31) 


As it has already been pointed out in Ref. 


53l |. looking at the correlation between and 


(/ / f) is quite useful to distinguish the four types of Yukawa interactions. 

In Fig. [H we show the tree level predictions on the Ake-A.kd plane (left panels) and Ake-Aku 
plane (right panels) in the four types of Yukawa interactions, where Akx = nx — 1- The subscripts 
E, D and U respectively represent the flavour independent charged leptons, down-type quarks 
and up-type quarks. In this plot, we take |a:| = 0.2, 0.14 and 0.028, and the sign of x is set to be 
negative (positive) for upper (lower) panels. As it can be seen, the predictions for the four types 
of Yukawa interacitons appear in different quadrants of the Ake-Ake plane. Therefore, at least 
from the tree level result, we can discriminate the type of Yukawa interaction in the THDM by 
looking at the measured values of Ake and Ake- 

In Ref. [^, one-loop corrected Yukawa couplings have been calculated in the four types of 
Yukawa interactions in the THDM. It has been clarified that the predictions in the four types of 
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FIG. 1; Tree level predictions on the (left panel) and Ake-^kjj (right panel) plane in the four 

types of Yukawa interactions. The black, blue and red curves respectively show the case of |a;| = 0.20 
[sin(/3 — a) ~ 0.98], |a;| = 0.14 [sin(/3 — a) ~ 0.99] and Jx] = 0.028 [sin(/3 — a) ~ 0.996]. The sign of x is 
taken to be negative in the upper hgures and positive in the lower figures. 

Yukawa interactions are well separated on the Ake-^kd plane at the one-loop level even if we 
scan the inner parameters under the constraints from perturbative unitarity and vacuum stability. 

B. Vacuum stability and perturbative unitarity 

A set of quartic coupling constants in the Higgs potential A 1 -A 5 is constrained by taking into 
account vacuum stability and perturbative unitarity as follows. 
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First, we require that the Higgs potential is bounded from below in any direction with a large 
ar held value. The sufficient condition to keep such a stability of the vacuum is given by [19l . 


sea 

lar 

00 

00 


Ai >0, A 2 > 0, \/Ai A 2 + A 3 + MIN(0, A 4 + A 5 , A 4 — A 5 ) > 0. 

Second, the perturbative unitarity bound 


(32) 


89l-l92l| is given by requiring that all the independent 


eigenvalues of the T matrix {i = 1-6) for the S-wave amplitude of the elastic scatterings of 
2 -body boson states are satished as 

1 


i»“±i < 5. 


(33) 


where each of is given by 




ai,± = 


1 


327r 

1 


= 


327r 

1 


3(Ai -|- A 2 ) ± y^9(Ai — A 2 )^ -|- 4(2A3 -|- A4)^ 
(Ai -|- A 2 ) ± (Ai — A 2 )^ -|- 4 A 4 

(Ai -I- A 2 ) ± (Ai — A 2 )^ -I- 4A| 


®4± = + 2A4 ± 3A5), 

’ ibvr 

05,± = -^4), 

« 6 ,± = 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


In Fig. [21 we show the allowed parameter region on the m 4 )-sin (/3 — a) plane (m<j> = mjj± = 
mA = rnu) from the constraints of vacuum stability and unitarity. It is seen that a large mass of 
additional Higgs bosons is allowed in a case with sin(,0 — a) ~ 1 . As another view of this hgure, 
we can extract the scale of the mass of the second Higgs boson from the precise measurement of 
Ky using Eq. (1271) . For example, if 1% deviation in the hVV coupling is found at future collider 
experiments, then the second Higgs boson should exist below about 800 GeV. 


C. The oblique parameters 


The S, T and U parameters proposed by Peskin and Takeuchi 


93 1 are modified in the THDM 


from those predicted in the SM due to the additional Higg s boson loop contributions and modified 
values of the SM-like Higgs boson coupling constants [60(]. We define the differences of S, T and 
U parameters as AS = S'thdm — S'sM, AT = Tthdm — Tsm and AU = C/thdm — Usm- These are 
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FIG. 2: The upper limit on the mass of additional Higgs bosons m$(= as a function of 

sin(/3 — a) for each hxed value of tan/3 in the case of cos(^ — a) < 0. The left regions from each curve are 
allowed by the constraints of vacuum stability and unitarity. 


calculated in terms of x defined in Eq. 


as 


AS = 


dvr 


- ^lnm^± 


+ 


AT = 


1 


+ 


\mh mh^ 

F^{Q\mA,rnH±) + 


fniA mz 


rriH rnn 

mfj± mA 

ruH ’ rnn 


+0(A 

/ ^mz mz^i I 


(40) 


m|,FA mlF^ / 

\mH mu J \ ruh mhj \ 


mu rrih 


mw rnw 


+ 1X1^ F/\ 


mz mz 


+ 4mwGA 



mh \ 

-4m|GA(— 

^)1 

\ 

\mw 

mw ^ 

\mz 

mzF 

1 


+ 0(x3), (dl) 


AU = 


\ ( mA mn \ 1, 2 77,// 2 ^ 

->- - - F^imz]mA,mH) 

dvr I \mH± mH±'' 3 


+ X 

+ 


/n^ i^\ / mH± mw \ _ p, mw \ 

^\mH ’ mn-' ^ ^ mh ’ mh / ^ v mh ’ mh / mu ' mu ) 


.mh mh 


mh mh 


mu mu . 


G'd^,^)-G'd^,^) +OF), 

\mz mz^ \mw mw^ 


(42) 
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where F^{my] 1711 , 1712 ) = [F^{mY; mi, m 2 ) — ^ 5 ( 0 ; mi, m 2 )]/my. The loop functions are given by 


F^{p‘^,mi,m 2 ) = [ hx [( 23 : — l)(mi — m|)+p^(2x — 1)^] In 

Jo 




2 -21 , ^2 


1-xl 


2 ) 


Fa(xi,X 2) = -(Xi - X 2 ) + 2 

^ L X 

n ( \ 1 , 1 + 2 1 + ®2 , 2 

Gi(I'l. X,) = 2 >“ - 1(13^ I"^'1 - ^(13^ I"^■2' 

’2(xf-xi)(l-xfx|) xf(xf-3) 2 , a;|(x|-3) 


1 


F^(xi,X 2 ) = - 


3[ (1 — x2)2(l — X2)2 (1 — x2)3 


In Xi + "' 2^2 In X 2 


(1 - x|)3 


Ga(xi,X 2 ) = 2 


1 — x| + 2 xf In xf ^ 1 — xf + 2 x| In x 2 


(1-X2)3 


(1-X2)3 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

h “ '"-2 

which gives zero in the case of mi = m 2 . Therefore, it is seen that AT becomes zero when x = 0 
and mA = or x = 0 and mu = Baii± is taken. 


where 


As = —x(l — x)p2 + xm\ + (1 — x)m 2 . 


In the case of p2 = g, the T 5 function is expressed by 

T5(0;mi,m2) = ;^(m2 + m^) + 


2 m? mn , m 2 
■ In —, 


2 ' ^ ’ m? — m? mi 


D. Flavour Constraints 


The mass of can be constrained from various B physics processes, because contributions 
from the SM IT-boson mediation are replaced by H^. In most of the cases, the constraint from 
s'y process provides the most stringent lower limit on m^± j^ . 1^ . In Ref. [^, the 


the b 


branching ratio of i? —)• Xg'j has been calculated at the next-to-next-to-leading order in the Type-I 
and Type-II THDMs. A lower bound has been found to be mi{± > 380 GeV at 95% conhdence 
level (CL) in the Type-II THDM with tan /3 > 2. A stronger bound for mjj± is obtained for smaller 
values of tan (3. On the other hand, in the Type-I THDM, the bound from 6 —)• 57 is important 
in the case with low tan/3; e.g., mH± < 200 (800) GeV is excluded at 95% CL in the case of 
tan/3 = 2 (1). When we consider the case with tan/3 > 2.5, the bound on m^± is weaker than 


the lower bound from the direct search at LEP, namely, about 80 GeV 


94l |. The similar bounds 


as those given in the Type-II and Type-I THDMs can be obtained in the Type-Y and Type-X 
THDMs, respectively, because of the same structure of quark Yukawa interactions. 

For a large tan /3 case, bounds from H —)• rz/ 


magnetic moment 


98 


95 


96| | , T —)• pi/u [96j, [971] and the muon anomalous 


99(1 can be more important as compared to the bound from 6 —)• 57 in the 
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Type-II THDM. For example, the lower limit on to be about 400 GeV is given at 95% CL in 
the case of tan/3 > 50 in the Type-II THDM 96]. 


For a small tan/3 case, the B^-B^ mixing is getting important to obtain a severe constraint 
on mf^± in the THD Ms. In the case of tan/3 = 1, m}j± < 500 GeV is exluded at 95% CL in all 


the types of THDMs 


lOOt j. This gives the stronger (weaker) bound than that from 6 —)• sy in the 


Type-II and Type-Y (Type-I and Type-X) THDMs. 


E. Direct searches for additional Higgs bosons at the LHC (7-8 TeV) 

The neutral Higgs bosons in the MSSM have been searched in the decay mode in the gluon 
fusion and bottom quark associated productions [^,38] using data with 7 TeV and 8 TeV of the 
collision energy and 4.9 fb“^ and 19.7 fb“^ of the integrated luminosity, respectively. Because the 
production cross section of the CP-odd Higgs boson from the bottom quark associated production is 
proportional to tan^ /3, high-tan/3 regions can be excluded by this process. For example, tan/3 > 10 
and tan /3 > 40 have been excluded at 95% CL for the fixed value of the mass of the CP-odd Higgs 
boson to be 300 GeV and 800 GeV, respectively 2. We can obtain a similar bound on tan/3 for 
a fixed value of tua in the Type-II THDM, because the structure of the Yukawa interaction is the 
same as that in the MSSM. Although the Hff coupling constant can be different in the Type-II 
THDM and the MSSM, we can achieve a similar value by taking sin(/3 — a) ~ 1, especially for the 
case with a rather large mass of the CP-odd Higgs boson in the MSSM. 

When sin(/3 —Of) / 1 is given, H —)• W^W~ jZZ decays can open in addition to the decay modes 
into a fermion pair. The search for the H —)• WW —)• ei'ni' signal has been performed 39|] in the 
range of 135 GeV < mn < 300 GeV using data with 8 TeV of the collision energy and 13 fb“^ 
of the integrated luminosity. The bound is presented in the mn-cosa plane for each fixed value 
of tan/3 in the Type-I and Type-II THDMs. In the Type-I THDM with tan/3 > 1, the strongest 
lower limit on mn is given to be about 220 GeV at 95% CL. On the other hand, in the Type-II 
THDM, similar bounds have been given as in the Type-I THDM. However, for a case with large 
tan /3, the excluded regions are shrinked due to an enhancement of fermonic decay modes such as 
H bh. 

In Ref. [^, H ^ hh and A ^ Zh decays have been searched in the THDMs with data of 
the collision energy to be 8 TeV and the integrated luminosity to be 19.5 fb“^. Multi-lepton and 
di-photon final states have been used for this search. The upper limit on the cross section times 
branching ratio has been presented for each of the processes gg ^ H ^ hh and gg ^ A ^ Zh ] 
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e.g., the upper limit of 8 (4) pb is given for the case of mn = 260 (360) GeV in the H ^ hh decay, 
while that of 1.6 (1.0) pb is given for the case of niH = 260 (360) GeV in the H —)• Zh decay. 
These bounds can be translated into the excluded regions on the cos(/3 — a)-tan/3 plane for given 
values of mu depending on the type of Yukawa interaction. 


F. Measurements of the Higgs boson couplings at LHC (7-8 TeV), and future collider 
experiments 


Both the ATLAS and CMS Collaborations have provided scaling factors for the Higgs boson 
couplings extracted from combined data of Higgs boson searches with y/s = 7 and 8 TeV and 


25 fb ^ of the integrated luminosity 




71]. Under assumptions of the universal scaling factors for 


fermions and vector bosons; i.e., Kp = Kt = Hj, = Kr and Hy = Kyr = k^, current data gives 


Ky = 1.15 ± 0.08, Kp = 0.99+[];?^, ATLAS [||, 


Ky = 1.01 ± 0.07, Kp = 0.87l[j;j3, CMS 




(50) 

(51) 


from the two parameters {Kp and Ky) fit analysis based on Ref. lOll j. The scaling factors for the 
loop induced Higgs boson couplings Kg and Ky have also been measured under the assumptions of 
Kp = Ky = 1, 


K, = 0.89;0:}1, 


Ky — l.iy_QJ2) 
K — 114 + 0 - 


ATLAS 11, 

CMS y. 


(52) 

(53) 


from the two parameters (k„ and Ky) fit analysis based on Ref. 101|. We can see that all the 


SM predictions {k^ = 1) are included within the 2-cr uncertainty of the measured scaling factors, 
where the current l-cr uncertainties of the scaling factors are typically of 0(10%). 

These scaling factors are expected to be measured more precisely at future collider experiments 
such as the HL-LHC and the ILC. In TABLE [III expected accuracies of the measurement for the 
scaling factors are listed at the LHC and at the ILC with several collision energies and integrated 
luminosities. 


HI. RENORMALIZATION 

We discuss the renormalization of the Higgs boson couplings, i.e, hZZ, hWW, hff and hhh 
at the one-loop level. In previous works, each part of the renormalized Higgs boson couplings has 
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Facility LHC 

HL-LHC 

ILC500 

ILC500-UP 

ILCIOOO 

ILClOOO-up 

yi (GeV) 14,000 

14,000 

250/500 

250/500 

250/500/1000 

250/500/1000 

/ Cdt (fb-i) 300/expt 3000/expt 250-b500 1150-bl600 250-b500-bl000 1150-bl600-b2500 


K~f 

5-7% 

2-5% 

8.3% 

4.4% 

3.8% 

2.3% 

Kg 

6-8% 

3-5% 

2.0% 

1.1% 

1.1% 

0.67% 

KW 

4-6% 

2-5% 

0.39% 

0.21% 

0.21% 

0.2% 

KZ 

4-6% 

2-4% 

0.49% 

0.24% 

0.50% 

0.3% 

Ke 

6-8% 

2-5% 

1.9% 

0.98% 

1.3% 

0.72% 

Kd = Kb 

10 - 13% 

4-7% 

0.93% 

0.60% 

0.51% 

0.4% 

Ku = Kt 

14 - 15% 

7 - 10% 

2.5% 

1.3% 

1.3% 

0.9% 


TABLE II: Expected precision on the Higgs boson cou pling s and total width at the 1-cr level from a con¬ 


strained 7-parameter fit quoted from Table 1-20 in Ref. 


102 |. 


been calculated. The one-loop corrected hZZ and hhh couplings have been evaluated in Ref. [7f 
in the Type-II THDM, and the hf] couplings have been calculated in Ref. Q in the four types 
of THDMs. 

We perform renormalization calculations based on the on-shell scheme which has been applied 
in Ref. However, it has been pointed out that there remains gauge dependence in the 


determination of the counter term of (3 in Ref. 


We thus construct a new renormalization 


scheme for f3 to get rid of the gauge dependence. As pointed out later in the paper, the gauge 
dependence is not completely removed, but shifted to a sector which does not contribute to the 
investigated couplings. 

First, we prepare a set of independent counter terms by shifting all the relevant bare parameters 
in the Lagrangian. We then give the renormalized one- and two-point functions which are written 
in terms of the contributions from IPI diagrams and counter terms. After that, we set the same 
number of renormalization conditions as the number of independent counter terms to determine 
them. 


For the determination of the counter term for the minimal subtraction scheme has been applied. 
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A. Parameter shift and renormalized functions 


We first perform the parameter shift of the electroweak sector and Yukawa sector as the following 


rriy —)■ rriy + 6my {V — W, Z), CJem Ckem + 
ruf ruf + 6mf, T \^2 —^ <^^ 1 ^ 2 , 

—y TP^ -j- oi —^ oi 6cx^ /? — y 13 -\~ 3^3^ — y A 

where ip = H^, A, H and h. The wave functions for the SM gauge bosons and 1T“ and the 
SM left (right) handed fermions tpi {'ipn) are shifted as 

^ (1 + w; ^ (^1 + ^6Zw^ w^, i^L/R ^ (1 + i^L/R- (55) 

We can then write down the renormalized two point functions for each particle. In the following, 
Bxy{p^) and n^y(p^) respectively denote the renormalized two point functions and the IPI dia¬ 
gram contributions for fields X and Y with the external momentum p^. The analytic formulae for 
the IPI diagram contributions are given in Appendix C. For the gauge boson two point functions 
IF+VF-, ZZ, 77 and the Z-'y mixing, we have 


^ww{p^) — n|y^(p^) — 5rn^ + 5Zyr{p^ — rn^), 

fizzip^) = n^^(p^) - (5m| + 6Zzip‘^ - m|), 
n^^(p2) = niPi(p2)+p2^z^, 

flz^ip"^) = nf^\p‘^) - SZz^ (^p‘^ - 

where 



bZz^ = Cy/Sy/(bZy/ — 5 Zb) = 2 ^ ^2 i^Zz ~ 3Z^). 

Cw - 


The renormalized fermion two point function is expressed by the following two parts: 


= Uffyip"^) + Uffyip"^), 


where 


Uffy{p^)=^ Ufj^y(p^)+6Z^ 


+ mf 


Uff^siP^) - 6Z^ 


n//A(p") =-7^5 uy-U(p') + ^zj, 


6mf 

ruf 


(56) 

(57) 

(58) 

(59) 


(60) 


(61) 


(62) 
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with 




(63) 


2 ’ ^ 2 

In Eq. ([62]), Hj-j y, and are the vector, axial vector and scalar parts of the IPI diagram 

contributions at the one-loop level, respectively. 

For the scalar sector, we first define shifts in the weak eigenbasis of the scalar fields: 



hi 

M, 



^odd 


W. 


. Wr, 


Z± 


w. 




(64) 


where Ze^em -^odd aiid Z± are arbitrary real 2x2 matrices. We then express shifts of the scalar 
fields in the mass eigenbasis as 


H\ H' 

I R{ — Soi)Zeven I 

h \h 


G^\ G^\ G^ 


R{-6I3)Z± 




(65) 


where we introduce Z^ven = R{—Oi)ZeveaR{a) and Zodd/± = R{~P)Zodd/±R{P)■ We define the 
matrix elements of them as follows: 


Zpvpjn — 


1 -I- ^Zh SCnh 
5GhH '^ + \Zh, 


Zodd — 


1 -I- 5Cga 
^Cag 1 + \Za . 


, = 


1 + \Zg± SCg+h- 
^Ch+g- l + 


( 66 ) 


We note that in Ref. [75[, the above matrices are chosen to be a symmetric form; i.e., SGffh = SGhu, 
SGga = ^Cag and 5Gg+h- = <5C'h+g-- la this paper, we do not take the symmetric form, and 
we use the additional degrees of freedom to remove the gauge dependence in the renormalization 
of 5/3 as it will be discussed in Sec. III-D. Finally, we can express the shifts of the scalar fields by 

1 -|- \6Zh SGffh + 

,SGhH — Sa l + ^6Zh 

1 -|- ^5.Zgo 5Gga + 5/3 
^SGag — 5/3 1 + ^6 Za 

/ 1 ^dZjj+ 6Gg+h- + 5/3 . . ^ 
y5C'jy+G- — 5/3 1 -b ^dZj^± 

For the scalar sector, we have the renormalized one-point function for h and H as 




fh = 6Th + Tj^\ fH = STh + I, 


( 68 ) 


18 




where 



The renormalized two-point functions are expressed as 

flhhip'^) = + [(P^ - ml)6Zh - 6ml] > (70) 

^hh{p‘^) = ^^fhip"^) + [(P^ - mjj)5ZH - 6mjj] , (71) 

nyiA(p^) = + [(P^ - m\)6ZA - 6m\] , (72) 

^H+H-ip"^) = ^H+H-ip"^) + [(P^ - rn]j±)6ZH± - 6m]j±\ , (73) 


and those of the scalar mixings are given by 


^Hh{p^) — n^(p^) + p^{6ChH + 6CHh) + ml{6a — 6ChH) ~ m\[{6a -|- SCuh)-, (74) 

nAG(p') = nAG(p') + p\6Cag + 6Cga) + m\{5(3 - 5 Cag), (75) 

nH+G-(P^) =n}f+g,_(p^) + p^{6Ch+g- +6Cg+h-) + mlj±{6fi — SCjj+q-), (76) 


where 


nir(p^) 

rlPI tJ2\ 


nff (/) + 


CpV 


SjjV 


nirA(p^) = n}rA(p^) + 

rlPI^„2 


IPI t‘2\ . Cq,(5Ti ^ Sa5T2 


+ 


c/^v spy 
slSTi cl6T2 

nyAp") = nS(p") + ^ + ' 


Cpv 


spy 


n 


H+H 


_(/) = n}P,V(p^) + 


sl5Ti cl6T2 


C0y 


+ 




spy 


nS(p^) = nS(p^) - (ih _ ih) , 

Vc/ju s/juy 

nAG(P^) = nAG(P^) + ^ [sin(/3 -a)TH- cos(/3 - a)r/,], 
n]|+G- (p^) = njf+G- (p^) + “ [®“(/^ -a)TH - cos(/3 - a)r/,]. 


(77) 

(78) 

(79) 

(80) 
(81) 
(82) 
(83) 


B. Renormalization conditions in the electroweak gauge sector 

The renormalization of the electroweak parameters can be done in the same way as in the SM, 
because the number of parameters to describe the electroweak observables are the same in the 
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THDM. This nature is also applied to models based on the SU{2)l x U{1)y gauge symmetry with 
p = 1 at the tree leveP. 

We apply the electroweak on-shell scheme based on Ref. 
counter terms in the electroweak sector; i.e., (Jaem, 5Zw and 5Zb- Therefore, we need 

the following five renormalization conditions to determine them: 


103l | to our model. There are five 


'Refiww {'^w) — 0) Ren2z(m|) = 0, 

An77(p^) 

dp^ ' p2=o 

= 0 , jji= j/2= me) = iejf,, 


= 0 , = 0 , 


(84) 

(85) 

( 86 ) 


where is the renormalized photon-electron-positron vertex. From the above conditions, we 
obtain 


dm^ — RenJy^(m(j/), dm^ — (j^z) ^ 


iPi/ 


5ae 




= n;ri(o)' - 

” cry m| 


5Z^ = -niPi(o)', = -^n^^^^(o) + 


2 —lPI/n^ I 


m 


where niP^(O)' = 


p 2=0 


swcw 

The other counter terms are also determined by 
rlPI/ 


dZz = ( 0 )' - ~ ~ 


civsw 

riPR 


mt 


-w 


9 ’ 


sz^v = -n‘Pi(o)' - ?£2llhhll + 




’ 




'-W 

L 


n^^(m|) n^L(mL) 


mr 


m 


w 


The counter term for the VEV Sv is also obtained through the tree level relation: 


^2 ^ '^w^w 


(87) 

( 88 ) 


(89) 

(90) 

(91) 


(92) 


as 


5v 

V 


Sw ~ ^w^wwi'^^w) I Cyp TTlPIm^/ I ^^7 


+ 


’VP 


m 


w 


’VP 


mt 


-n;(;^( 0 )' + 


cvp mt 


(93) 


® When we discuss models without p — 1 at the tree level such as models with isosipin triplet scalar fields, one 
additional input parameter is required to express the electroweak sector. Therefore, we need an additional renor¬ 
malization condition to determine the extra counter term associated with the parameter. In the model with a 
y = 0 Higgs triplet field, the renormalization of electroweak parameters has been discussed in Refs. [^ . 1^ . 
Furthermore, in the model with a T = 1 Higgs triplet field, that has also been discussed in Refs. (b^. F^. 
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We here note that the fermion-loop contribution to is given by 


n77'(o)' = E "^/)’ (94) 

where Qf is the electric charge of a fermion /, A^/ is the color factor: nI = 3 (1) for / being quarks 
(leptons), and A is the divergent part of the loop integral as defined in Eq. (IB23jl in Appendix B. 
In order to avoid to input the light quark masses, we can use the following relation obtained from 
Eqs. (I58l) and (f88l) 


rur. 


nii;i(o)' = — nipi(m|) - n,,(m|) 


= Aoe 


m% 


77 


(95) 


where Actem is the shift of the structure constant that we can quote the experimental value. In 
the right hand side of the above equation, the light fermion mass dependence in is 

of order so that we can neglect it. 


C. Renormalization conditions in the Yukawa sector 

In the Yukawa sector, there are three counter terms 5mf, 6Zy and To determine them, 

we impose the following three conditions for the fermion two point functions H: 


n//y(?^/) = o, ^n//,y(p 2 ) 


= 0 , 4-AffAip‘^) 

p^=my d'^ ’ 


p^=m?^ 


= 0 , 


(96) 


we obtain 


5m 


f - rrlPI 




mt 


5Zy — —(ttt-j) — 2mj 




O 9 

j 


d 


5Z{ = + 2m^^n}Pi^(p2) 


p^=m} 


(97) 


D. Renormalization conditions in the Higgs potential 


There are totally 21 counter terms in the Higgs potential, namely, the counter terms for two 
tadpoles and STh, four mass parameters {ip = A, H and h), two mixing angles 6 a 
and 5(3, four wave function factors 6Z^, six wave function mixing factors 5Cij, and 6M‘^ First, 

^ In addition to them, there are two more counter terms 3Zq± and SZqo. However, they do not enter the following 
discussion. 
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we impose two tadpole conditions at the one-loop level, i.e., 


Th = TH = 0 . 


(98) 


We then obtain 


6n = 5 Th = -r 


H 


Second, eight on-shell conditions for the two-point functions: 

= 0 , 

= 0, for ip = A, H and h, 

which determine the following eight counter terms 

4<5ri , cl5T2 




C/3V 




S/3V 


,2 _ rrlPI ^ S^5T2 


6 m H — + 


6 m\ = + 


Cpv S0V 

sl5Ti cl6T2 


c^v 


+ 


■'/S'- 


spy 


Sm'jp± = n]P^_(m^±) + 


slSTi ^ 46 T 2 


c/sv 


S/3V 


and 


d 




(99) 

( 100 ) 

( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 

(106) 


'P'P 

Three counter terms 6 a, 6ChH and dCuh related to the mixing between the CP-even scalar 
states are determined by imposing the following three conditions 


n_ffh(w.|) = = 0, 6 ChH = dCuh = 6 Ch- 


They give 


6 a = 


1 


6Ch = 


1 




6T\ _ 

CpV S0v)\ ’ 


(107) 

(108) 
(109) 


2 {m]j - ml) 

Three counter terms 6(6, 6Cag and 6Cga related to the mixing between the CP-odd scalar states 
are determined by three conditions. Similar to the CP-even sector, we first impose the following 
two conditions as 


nAG(o) = n^G(”ii) = 0- 


( 110 ) 
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We then obtain 


5/3-5C'^G = -^n3iG(0), 5l3 + 5CGA = -\ti^Akm\). (HI) 

m\ m\ 

In order to determine three counter terms, we need to impose one more renormalization condition 
in addition to that given in Eq. (IllOp . This third condition can be used to remove the gauge 
dependence in h/I which was already mentioned in the beginning of this section. To define such a 
condition, we separate into the gauge dependent (G.D.) part and the gauge independent 

(G.I.) part as 




( 112 ) 


Then, we imposed the third condition as 


1 




(113) 


Using Eq. (Illip . the remaining two counter terms are also determined: 

1 


5Cag = - 


SCga = - 


2 m\ 


1 


2 m\ 


nAG("ii)lG.i.-2nAG(0)|G.D. 
n)4G("ii)|G.L + 2nAG("ii)|G.D. 


(114) 

(115) 


We note that in n(^^(0) only the G.D. part is survived; i.e., n35^(0) = n3j^(0)|gjj . As it can 
be seen in Eqs. (jll4p and (jll5p . there still remains the gauge dependence in 5Cag and SCga- 
However, they do not appear in the following calculations for the renormalization of the Higgs 
boson couplings. Instead of applying the above renormalization scheme for h/3, we can apply the 
MS scheme in which the gauge dependence can also be removed at the one-loop level as discueed 
|. In the following discussion, we apply the renormalized tan/I determined by Eq. (Ill.Sp . 


in Ref. 


The above A-G^ mixing can be replaced by the mixing between A and the physical Z boson 
by the help of the Ward-Takahashi identity; i.e., the condition IlyiG(^A) = 0 is equivalent to that 
of vanishing renormalized A-Z mixing; i.e., tizAi^nn?^ = 0, which can be defined in the following 
way. The Z-A mixing is obtained from the kinetic term: 


Ain = mz{d^,G^)Z^^ + • • • ^ mz{5p + 5CGA){d^A^)Z^^ + • • • 


(116) 


The renormalized Z-A mixing H^^ = —ip^IizA{p^) is then expressed by 


IizA{p^) = mz{5p + 5Gga) + 


(117) 
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SQ set set 

Type-1 

— ^ [cot (36 P + tan aSa) — ^ {cot (36(3 + tan aS a) — ^ {cot (36(3 + tan a 6 a) 

Type-Il 

— ^ {cot (36(3-\-tan a 6 a) — ^ {tan (36(3 + cot a 6 a) — ^ {tan (36(3-\-cot a 6 a) 

Type-X 

— ^ {cot (36(3 -\- tan a 6 o() — ^ {cot (36(3 + tan a 6 a) — |^(tan/35/3-I-cota(5ck) 

Type-Y 

— ^ {cot (36(3 + tan a 6 a) — ^ {tan (36(3 + cot 060 ) — ^ {cot (36(3 + tan a 6 a) 


TABLE III: The counter term for the mixing factors in Yukawa interactions. 


where is the incoming momentum of A. The IPI diagram contribution to the Z-A mixing 
is given in Appendix. Because of the relation /mz, the con¬ 
dition = 0 can be replaced by n 2 ^(m^) = 0. Therefore, Eq. (|113l) is rewritten as 

1 


<5/3 = - 


2 mz 




(118) 


We note that the numerical difference between in our scheme and in the previous scheme applied 
in Ref. s is negligibly small as long as we discuss the case with sin(/3 — a) ~ 1 or x ^ 1. 

Two counter terms SC^+g- ^Cq+h- for the mixing between the singly-charged scalar 
states are determined by requiring the vanishment of the mixing between and at = 0 


and p^ = m 


H±- 


ni/+G-(0) — 'Ah+g-(j^‘h±) — 0- 


(119) 


We obtain 


1 


5Ch+g- ^n]:f.J^_(0), 5Cg+h- = -S/3 - ("^h±)- 


1 


ipi 


m 


W± 

Until here, we did not discuss the determination of (5M^. As adopted in Ref. [75|, we apply the 
minimal subtraction scheme for 5M^, where it is determined so as to absorb only the divergent 
part in the hhh vertex at the one-loop level, that is 

j. Sin zp 

( 121 ) 


( 120 ) 


IGvr^u^ 
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IV. ONE-LOOP CORRECTED HIGGS BOSON GOUPLINGS 


A. Analytic expressions 

In the previous section, all the counter terms are determined by the set of renormalization 
conditions. Now, we can evaluate the one-loop corrected Higgs boson couplings hWW, hZZ, hff 
and hhh. In addition to the above couplings, we also give formulae for the loop induced decay 
rates h —)■ 77 , h —)• Z 7 and h —>• gg. 

The renormalized hVV, hff and hhh vertices are expressed as 


2 2 2 \ 
1,P2,Q ) 

2 „2 „ 2 \ 


^iff{Pl,P2 



( 122 ) 

ffT+Kff+^^h}?iplpl9% 

(123) 

'hhh + ^^hhh + ^jKhiPl^P2^ 9^)) 

(124) 


where , dT hxx find T j^x ® contributions from the tree level, the counter terms and the 
IPI diagrams for the hXX vertices, respectively. In the above expressions, pi and p 2 {q = pi +P 2 ) 
are the incoming momenta of particle X (outgoing momentum for h). 

For the hVV and hff vertices, the indices i and j label the following form factors: 


_ -pl I 1^2 P 1 P 2 I ,-p 3 ^^upcrP^pP‘^0- 

^hVV - ^hvv9 + ^hVV—r 0 ’ 

ffLy Ifly 

= ^hff + Tsf f + f)2 

+ +i>2ib^hjf + + i>ii>2ib^hjf- 


(125) 


(126) 


The tree-level contributions are given as 


-\l,tree 

hVV 


• m ^ 
- — sm(p — a), 

V 

3,ti 
hVV 


ptree _ fffJLcf 

'■hff - 


p2,tree _ p3,tree _ pj,tree _ p, ( ■ / c\ 

hVV — ^ hVV —^hff “ U T 


hff 


ptree _ 

^ hhh 


- 6 A 


hhhi 


(127) 


The counter-term contributions are 

sin(/3 — a) 


1 _ 2ml r 

hVV — - 


V 


dniy 


1 


5v\ 


— 9 — + SZy + —6Zh -J + cos(/3 — oi){6f5 + 5Ch') 

rriy 2 V J 


AF'^ — cf 

hff - ^ Zh 


f f 

6mf 6v in ^ X 

-h 6Zy + —SZfi H- J- H- j{5Ch + Sa) 

"V ” 2 (I (I 


^^hhh = 6 


S^hhh + 2^^h + ^Hhhi^Oi + 6Ch) 


^^hVV — ^^hVV — ^^hff — U 7 ^ S), 


hVV 


hff 


(128) 


25 














where 


6v 1 

o^hhh — —^hhh -1- 


V V sin 2/3 


cos^(;8 — a) cos(a + j3)5M‘^ 


1 


Av sin 2/3 


[cos(3a — /3) + 3cos(a + /3)] 5m\ 


1 


+ TT ~ “) 
2v 


sin 2a, , , o, , ,n 

sin 2(3 


{mi - M^) + 


da 


+ 


1 


4r; sin^ 2/3 


cos(/3 — a) [(4 + 4 cos 2a cos 2/3 — 2 sin 2a sin 2l3)mi 


— (5 — cos 4/3 + 4 cos 2a cos 213 — 2 sin 2a sin 2/3)M^](5/3. (129) 

The counter terms 5^^ appearing in the Yukawa couplings are expressed in terms of <5/3 and 6a as 
listed in Table m We dehne the renormalized scaling factors in the following way: 
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The momentum q'^ is fixed to be {my + m/i)^, and {2mh)‘^ for ky, k^ and respectively, in 
the following discussion. 

The deviations in the renormalized Higgs boson couplings are approximately expressed by keep¬ 
ing the non-decoupling effects of extra Higgs bosons and top and bottom masses dependence 
{mA — mu is assumed) as 
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where = 2 (1) for {H, A). We can see that there appears the term 

rr?^lv^ (l — M'^jrn?^ in A/ty which comes from the counter term bZh\ i.e., the derivative of 
the h two point function given in Eq. (jlOOp . When we consider the case with this term 

gives the quadratic power like dependence of the mass of additional Higgs bosons. This corresponds 
to the case where the masses of the additional Higgs bosons, which is expressed schematically as 
m\ = XiV^ + M^, mostly come from the Higgs VEV v. In such a situation, it is known that the 
decoupling theorem does not work. On the other hand, if we consider the case of ;:$> the 
amount of /S.kf is reduced as l/m| according to the decoupling theorem. The same contribution 
from 6Zh is also seen in Akj (/ = T,c,b,t) through the term Aky- Notice here that there are 
additional terms proportional to the top or bottom quark masses in Akb and Akt- Apart from Aky 
and Akf, let us discuss the expression of Akh- There appears the term rn^/ (l — M'^ 
which comes from the additional Higgs boson loop contributions to the IPI hhh diagrams. When 
we consider the non-decoupling case; i.e., it gives the quartic power like dependence of 

771$. Similar to the case in Any, this effect is decoupled by l/m| when is taken. 

Similarly, the decay rates of /i —)• 77 and h ^ gg are expressed in terms of x (x <C 1) as 
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r(h ^ gg) ~ 




^(1 +^qX- ^)If 


(138) 


where Ip and lyr are the loop functions. The exact expressions for the decay rates for h —>• 77 , 
h —)• Z 7 and h ^ gg are given in Eqs. ()C58I) . (IC59D and (jCGOp in Appendix C, respectively. In 


Eq. ()137p . the first term in r(h —)• 77 ) proportional to (1 — M^/m|^±) is the charged Higgs boson 
loop contribution. When we take the limit of —)• 0, this term approaches to the constant —1/3. 

This can also be understood as the consequence of the non-decoupling effect of the charged Higgs 
boson loop contribution, but it is not like the quartic (quadratic) power like dependence as seen in 
Akh {Aky and Akf). 
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FIG. 3; Deviations in the scaling factors for hVV (upper left), hbb (upper right), h'-fj/hZ^ (bottom left) and 
hhh (bottom right) at the one-loop level as a function of m^{= 'm^± = = m^) in the case of sin(/3—a) = 

1 and tan/3 = 1 The black, blue and red curves respectively show the cases of VXv‘^{= — IvP) = 150, 

300 and 400 GeV. 

B. Numerical evaluations 


In the following, we show numerical results for the Higgs boson couplings at the one-loop level. 
We use the following inputs 94l |: 


mz = 91.1875 GeV, Gp = 1.16639 x 10"® GeV"^ = 137.035989, Aaem = 0.06635, 

rrit = 173.07 GeV, mj, = 4.66 GeV, rric = 1.275 GeV, m-r = 1.77684 GeV, 

ruh = 126 GeV. (139) 


We first show the case of the SM-like limit a; = 0. In this case, the deviations in the Higgs 
boson couplings purely comes from the additional Higgs boson loop effects. We note that the 
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FIG. 4: Deviations in the scaling factors for hVV (upper left), hff (upper right), hjj/hZ^ (bottom left) 
and hhh (bottom right) at the one-loop level as a function of m^{= = rrijj) in the case of 

= 0, sin(/3 — a) = 1 and tan /3 = 1. 


tan/? dependence in the renormalized scaling factors appears only in kf. We take all the masses 
of additional Higgs bosons to be the same; i.e., m}j± = ruA = ruH (= ?7i$) for simplicity. 

In Fig. [3l we show the decoupling behavior of additional Higgs boson loop contributions to 
the Higgs boson couplings. The upper-left, upper-right, lower-left and lower-right panels re¬ 
spectively show Aky, Akh, and Akh as a function of for several fixed values of 

VXv^ (= m\ — ) in the case of tan (5 = 1. We can see that all the deviations approach 


to zero in the large mass region due to the decoupling theorem 


28|. 


In Fig. m we show the deviation in the Higgs boson couplings Aky (upper-left), Akef (upper- 
right), Ak^!^^ (lower-left) and Akh (lower-right) as a function of m^. We take = 0 and 
tan /3 = 1 for all panels. In this case, the magnitude of deviations increase when becomes 
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Set A 

Set B 

Set C 

Set D 

Set E 

Akv 

-2% 

-2% 

-2% 

-1% 

-0.4% 

AKt 

-bl8% 

-bl0% 

+5% 

-bl8% 

+18% 

Ak}, 

-bl8% 

-bl0% 

+5% 

+ 18% 

+18% 


TABLE IV: Benchmark sets for the central values of measured scaling factors for the hVV, hbb and Htt 
couplings. The expected 1-cr uncertainties for each scaling factor at the HL-LHC and the ILC 500 are shown 
in Eq. (I140j) . 


larger due to the non-decouipling effect of the extra Higgs boson loops except for Ak 


2 

7 /^ 7 ' 


V. DETERMINATION OF INNER PARAMETERS FROM THE HIGGS BOSON COU¬ 
PLING MEASUREMENTS 


In this section, we investigate how we can fingerprint the THDMs using the one-loop corrected 
Higgs boson couplings and also future precision measurements of these couplings at the HL-LHC 
and the ILC. We carefully see how the tree level analysis for the model discrimination discussed 
in Sec. II or in Ref. 53| can be improved by the analysis with radiative corrections. Furthermore, 


we demonstrate how the inner parameters such as x, tan (3 and masses of additional Higgs bosons 
can be extracted from the measurement of the couplings for the Higgs boson h. In our analysis 
below, we assume that the deviations in scale factors of the Higgs boson couplings are measured 
as expected in Table IIVI We also assume that the SM values of these coupling constants are well 
predicted without large uncertainties which mainly come from QCD corrections®. 

Let us suppose that Any, Ak-t and are measured at the HL-LHC and the ILC500. We 
consider five benchmark sets for the central values of (Any, Ak^-, Aki,) as listed in Table HVl Set A 
is the typical case where Yukawa couplings deviate from the SM values rather significantly (18%) 
with a relatively large deviation in the hVV couplings (—2%). Set B and Set C correspond to the 
cases with smaller deviations in Yukawa couplings with the same deviation in gauge couplings as 
Set A. Set D and Set E do to the cases with smaller deviations in gauge couplings with fixing the 


® According to Refs. [lOdl . Iiosl ]. the current uncertainty of the bottom Yukawa coupling hbb due to the QCD correc¬ 
tions i s 0.7 7% in the SM. This uncertainty could be reduced in future studies using the lattice calculation up to 
0.10% [l05| which is better than the expected accuracy of the measurement of the hbb coupling at the ILClOOO-up 
as listed in Table [IT] (0.4%). 
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same deviation in Yukawa couplings as Set A. According to Table HIl the l-a uncertainty for these 
scaling factors are given as 

la(Kv),o-(Kf,),<T(K-r)] = [2%, 4%, 2%], for HL-LHC, 

[(t(kv),(T(K fe),( t(kt-)] = [0.4%, 0.9%, 1.9%], for ILC500. (140) 

From the tree level analysis in Fig. [H these benchmark sets indicate that the Higgs sector is the 
THDM with the Type-II (Type-I) Yukawa interaction assuming x ~ cos(/3 — a) < 0 (x > 0). In 
order to further discriminate Type-I or Type-II, we need additional information to determine the 
sign of X such as the measurement of Akc, namely, if Akc is given to be a negative (positive) value, 
then we can completely determine the Yukawa interaction to be Type-II (Type-I). In the following, 
we consider the case of Akc < 0, so that we assume the case of the Type-II THDM. 

For all Set A to Set E, we survey parameter regions in which values of k’s are predicted around 
the central values within the l-cr uncertainty expressed in Eq. (11401) by scanning the inner pa¬ 
rameters X, tan /3, (= rnjj± = rrij^ = m^) and in the Type-II THDM. We also take into 

account the constraints from vacuum stability and perturbative unitarity in order to constrain the 
parameter space. The scanned regions for tan (3 and are taken as tan /3 > 1 and > 300 
GeV, respectively. Values of the other parameters and x are scanned over ranges which are 
enough wide to obtain the maximally allowed parameter spaces. 

In Eig. [5l we show the allowed parameter regions on the x-tan/3, x-m^, rn^-C, and m^-tan/3 
planes from the left to right panels, where we define 

?fi$ = m$C- (141) 

The parameters x and ffi^ give deviations of the Higgs boson couplings by the mixing effect and 
the loop effect, respectively. Notice that the scale of rh^ corresponds to the mass of the extra Higgs 
boson when = 0. The physics meaning of C is to measure the magnitude of non-decouplingness 
of the loop effects of extra Higgs bosons. If ^ is unity, we have = 0, while if C < 1 with 
nonzero value of (> 0), the mass of the extra Higgs bosons partially comes from so that 
the non-decouplingness is smaller. The central values of Ak’s are chosen from Set A, B, C, D and 
E from the upper to bottom panels. The blue and red points correspond to the region within the 
l-fj uncertainty at the HL-LHC and ILC500, respectively, from the central value in Table HVl 
For Set A in Fig. [5l let us first explain the behavior of the red points on the x-tan (3 plane. In 
this case, —2.4% < AKy < —1.6% is allowed at the ILC500, which can be explained by taking 
-0.22 < X < —0.18 at the tree level from the expression of Any ~ — x^/2. At the same time, 
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FIG. 5: Scatter plots for Set A, B, C, D and E from upper to bottom panels. The cyan and red points satisfy 
the benchmark sets within the l-cr uncertainty at the HL-LHC and ILC500 given in Eq. (I140|) . respectively. 
Eor the panels shown in the second and the third columns, the vertical axis m<i> and C are respectively 
defined by = m$(l — and C = 1 — 


both Akt and AKb are approximately given by —xtan/3 in the Type-II THDM at the tree level, 
so that tan/3 is determined by a fixed value of x from tan/3 ~ —AK^jblx, which is around unity 
if we take the central value of Akv and AK^/b- In fact, by looking at the top-left panel in Fig. [U 
the above mentioned values of x and tan/3 are allowed. However, the actual allowed region of x 
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inducing radiative corrections is about from —0.22 to —0.12 which is wider than the allowed region 
estimated at the tree level. This can be understood by taking into account the additional Higgs 
boson loop contributions to Ky at the one-loop level. The approximate formula for Aky is given in 
Eq. (11311) . where the second term in the right hand side corresponds to the one-loop contribution. 
The point here is that the sign of one-loop effect is negative, and it is proportional to the factor 
Therefore, the allowed region above x ~ —0.18 is explained from the one-loop contribution 
with a non-zero value of On the other hand, the one-loop correction to Kr is given by the same 
form as for Ky as given in Eq. (I132p . so that the difference Ak^ — Aky is approximately given 
by the same form — xtan/3 as that given at the tree level. Now from the measurement, since the 
difference is determined with the uncertainty, —x tan (5 is also fixed at the one-loop level. We thus 
can understand the shape of the allowed region of this plot. Although for Akb the top quark, the 
bottom quark and loop diagrams give an additional contribution as shown in Eq. (I134p . this 
is not so significant in the scanned regions. As a consequence for Set A, when the measurement 
at the ILC500 is assumed, the allowed value of x and tan can be determined to be about from 
—0.22 to —0.12 and from 1 to 2, respectively. On the other hand at the HL-LHC, Akv = 0 is 
included within the l-cr uncertainty. Thus, x ~ 0 is still allowed, so that the value of tan (3 is not 
determined at all because of the relation tan/3 ~ —AKT-jbjx. In addition, we can only extract the 
lower limit of x to be about —0.22. 

Next, we discuss the behavior of the second panel for Set A in Fig. [5l As we mentioned in the 
above, the vertical axis fh^ measures the size of one-loop contribution to the deviation in the Higgs 
boson couplings. At the ILC500, in the region with x ~ —0.20, the value of fh^ is determined to 
be a smaller value, but fh^ ~ 0 is not included because of the constraint from vacuum stability. 
This can be understood that the deviation from the tree level mixing is dominant in this case. 
On the other hand, when the value of x approaches to zero, a sizable value of is extracted, in 
which the deviation driven by the one-loop contribution becomes more important to compensate 
the reduced contribution from the tree level mixing. In addition, the upper limit of to be about 
450 GeV is determined by the constraint from perturbative unitarity. At the HL-LHC, although 
the blue plots are spread over the region with x ~ 0 as we observed in the x-tan/3 plot, the upper 
and lower limit of is given by the constraint from unitarity and vacuum stability, respectively. 

The third panel for Set A in Fig. [5] shows the allowed region on the m^-C, plane, where C, is the 
parameter indicating the non-decouplingness of the extra Higgs bosons. For Set A, the allowed 
regions for ILC500 are shown by the red points while those for HL-LHC by the blue points. There 
are upper and lower bounds for C, for each value of They are crossed at around = 850 GeV 
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which corresponds to the upper bound of the mass of extra Higgs boson. The region of C is from 0.2 
to 1.4 at = 300 GeV. The region of C > 1 corresponds to < 0, where non-decoupling effects 
are effectively large. The exclusion of C < 0.2 means that there must be some non-decoupling loop 
effects of extra Higgs bosons in order to explain this benchmark point. At the HL-LHC, the similar 
behavior can be observed. However, <^ = 0 is still allowed, so that we cannot say something about 
the non-decoupling effect. 

The last panel for Set A in Fig. [5] shows the allowed regions on the m,j-tan/3 plane. At the 
ILC500, tan (3 can be determined to be less than 2, and the upper bound of the mass of the extra 
Higgs bosons are obtained to be less 850 GeV, while at the HL-LHG, tan/3 is undetermined and 
only the upper bound of the mass of the extra Higgs bosons is obtained. 

The panels shown in the second and third rows in Fig. [5] display the allowed parameter regions 
for Set B and Set G, respectively, where the central value of Akt-{= Ak^) is taken to be smaller than 
that of Set A, while Any is taken to be the same. By looking at the panels for the x-tan /3 plane, 
we can see that a smaller value of |a:| is preferred as compared to the case for Set A. Furthermore, 
a smaller value of tan/3 is favored in addition to a smaller value of \x\ as seen in the result at 
the ILG500. These tendencies can be understood in such a way that the deviations in Yuakwa 
couplings are proportional to —xtan/3 at the tree level. Because of the smaller value of \x\, the 
deviation in Ky cannot be explained only from the tree level contribution, so that the one-loop 
effect is necessary to compensate the tree level contribution. That is the reason why the red points 
in the second and the third panels for Set B and Set G are given in the upper region which does not 
include fh^ ~ 0 and C — 0. Therefore, the non-decoupling effect can be extracted at the ILG500 
for these two benchmark sets. From the results of ILG500, the upper limit on is extracted to 
be about 950 GeV and 800 GeV for Set B and Set G, respectively. 

The panels shown in the fourth and fifth rows in Fig. [S] display the allowed parameter regions 
for Set D and Set E, respectively, where the central value of Any is taken to be smaller than that 
of Set A, while Anr (= Ak^) is taken to be the same. From the red points in the left panels, it is 
seen that the values of smaller |3:| and larger tan/3 are allowed, which can be explained by the tree 
level formulae of Any = —x^/2 and Ak^/ij = —xtan/3. For Set E unlike the other benchmark sets, 
values of x and tan /3 are not well determined even at the ILG500, because Any ~ 0 is included 
within the 1-a uncertainty of ILG500. The extraction for m^, and is done from the ILG500 
as 50 < < 300 GeV, 0.1 < C ^ 1-1 GeV and < 850 GeV for Set D and 0 < fh^ < 200 GeV, 

0 < C < 0.7 GeV and < 800 GeV for Set E. 

Up to now, we have discussed the extraction of the inner parameters from the three experimental 
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FIG. 6: Scatter plots for Set A with the additional constraint from = 0.98, 1.00 and 1.02 for upper, 
center and bottom panels. The 1 -ct uncertainty of Kj is assumed to be 2% as expected at the HL-LHC. The 
cyan and red points satisfy the benchmark sets within the 1-sigma uncertainty at the HL-LHC and ILC500 
given in Eq. (11401) . respectively. For the panels shown in the second and the third columns, the vertical axis 
and ( are respectively defined by = m$(l — A'P/m%) and C = 1 — M'^ /m%. 

inputs; i.e., A^ty, Ak,- and Ak;,. In Fig. (G] we show how the extraction can be improved by adding 
information of in addition to the above three inputs. The panels shown in the first row are the 
same as those shown in the first row in Fig. (5] which are displayed in order to compare the results 
with The panels displayed in the second, third and fourth rows respectively show the allowed 
region for Set A with the central value of of 0.98, 1.00 and 1.02 within the l-cr uncertainty of 
±2% as expected at the HL-LHC (see Table Hit) . Because the accuracy of the measurement of n-y 
at the ILC500 is not better than that of the best value at the HL-LHC, 2%, we also use 2% for 
the analysis at the ILC500. As we see Eq. (jl37p . the loop contribution to the decay rate of 
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the /i —)• 77 mode gives a different dependence of the non-decouplingness from that in Aky and 
Akf, which is not proportional to but proportional to C, so that the non-decouplingness C 
can be expected to be extracted more precisely depending on the measured value of In fact, 
we can observe that ( is determined more precisely to be 0.5 < C ^ 1-0) 0-25 < C ^ 1-1 and 
0-2 < C ^ 0.5 at the ILC500 for the cases with the central value of = 0.98, Kj = 1.00 and 

= 1.02, respectively, as compared to the case without (0.2 < C ^ 1.2). The determination 
of is also improved, because m<i> is given as a function of We note that smaller values of C 
and rh^ are favored in the case of the larger central value of k^, because the loop effect gives 
a destructive contribution to the W boson loop contribution. 

In Fig. [71 we also show the allowed parameter region with additional information of for Set 
D. Similar to the results in the previous figure, C and fh^ are well extracted as compared to the case 
without displayed in the first row in Fig. [71 For example, ( is determined to be 0.3 < C ^ 0-8, 
0-1 < C ^ 0-6 and 0.1 < C ^ 0.6 for the cases with the central value of = 0.98, = 1.00 and 

K-y = 1.02, respectively. 

VI. DISCUSSIONS AND CONCLUSIONS 

We have calculated radiative corrections to a full set of coupling constants for the Higgs boson 
h at the one-loop level in the THDMs with the four types of Yukawa interactions under the softly- 
broken discrete Z 2 symmetry. These couplings are evaluated in the on-shell scheme, in which the 
gauge dependence in the mixing parameter which appears in the previous calculation is consistently 
avoided. We have shown the details of our one-loop calculations, and have presented the complete 
set of the analytic formulae of the renormalized couplings. We then have numerically demonstrated 
how the inner parameters of the THDM can be extracted by the future precision measurements of 
these couplings at the HL-LHC and the ILC. 

We have found that the inner parameters of the THDM can be determined to a considerable 
extent as long as Ky will be measured with the deviation about 1%. The extraction of the inner 
parameters using the ILC500 is much better than that using the HL-LHC. That is mainly due to the 
good accuracy of the hVV coupling measurement at the ILC500 whose uncertainty is expected to 
be less than 1%. Although we have only demonstrated the results for Set A to Set E assuming the 
true Higgs sector is of the Type-H THDM, the similar analysis can be performed straightforwardly 
in the other types of THDM or the other extended Higgs sectors, and the extraction of inner 
parameters is expected to be attained as well in these models. Our study given in this paper 
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FIG. 7: Scatter plots for Set D with the additional constraint from K-y = 0.98, 1.00 and 1.02 for upper, 
center and bottom panels. The l-cr uncertainty of is assumed to be 2% as expected at the HL-LHC. The 
cyan and red points satisfy the benchmark sets within the 1-sigma uncertainty at the HL-LHC and ILC500 
given in Eq. (11401) . respectively. For the panels shown in the second and the third columns, the vertical axis 
771$ and C are respectively defined by 7fi$ = m$(l — M'^/m\) and C = 1 — M'^/m\. 

shows that the numerical evaluation of the Higgs boson couplings at the one-loop level in extended 
Higgs sectors is essentially important to indirectly determine the structure of the Higgs sector by 
using the future precision data. In addition, it also shows that in addition to the HL-LHC where 
especially /lyy can be measured precisely future lepton colliders such as the ILC are absolutely 
necessary for our purpose of determining the structure of the Higgs sector from the measurement 
of the coupling constants of the discovered Higgs boson h. 

Although we have discussed fingerprinting by using Ky, Kt, Hb and the information of Kc, 
and Kh is also important to determine the Higgs sector more deeply. In particular, the measurement 
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of the top Yukawa coupling is important not only to determine the nature of the top quark, the 
heaviest matter particle, but also to test the new physics scenarios based on the composite models. 
The measurement of the hhh coupling is essentially important not only to determine the nature of 
the Higgs potential but also to test, for instance, the new physics models with strongly hrst order 
phase transition. Although at the HL-LHC the cross section of the double Higgs production process 
is expected to be measured at a few times 10% it seems to be hopeless to extract the information 


of the hhh coupling sufficiently accurately. On the other 
hhh coupling can be measured with the 13% accuracy 59l . 


rand , at the ILC with y/s = 1 TeV the 


106l | , which is sufficient precision to test 


the strong first order phase transition which is required for successful electroweak baryogenesis. 

We conclude that the combination of the future data for all kinds of the couplings for the Higgs 
boson h and their theory predictions with radiative corrections in various extended Higgs sectors 
is a promissing way to determine the structure of the Higgs sector and further to access new 
physics beyond the SM, even if a new particle was not directly discovered in the future experiments. 
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Appendix A: Higgs boson couplings 

Erom the Higgs kinetic term, we obtain the two types of the trilinear couplings; i.e., Gauge- 
Gauge-Scalar, Gauge-Scalar-Scalar, and quartic Gauge-Gauge-Scalar-Scalar type couplings. These 
couplings can be expressed as 

^ = + Vi - 4’id^4'2)yti + g<i>i<i>2ViV29^''4 >i4>2Vi^V2u h— • (Al) 

The coefficients grj,ViV 2 ^ 9 <f>itf> 2 V 9(l>icl>2ViV2 listed in Table IVTl where we use g^ = g/cy^ in 
this table and below. Throughout Appendix, we use the shortened notation of the mixing angles, 
= sin(/3 — a) and c^_o = cos(/3 — a). 

Erom the Higgs potential, we obtain the scalar trilinear and the scalar quartic couplings. When 
we use the following notation for these couplings 
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2 

VCjS—a 

hZ^Zi, 

al 

HZ^Z, 

^VCp-a 

G^Z^W^ 

aaz „,„2 

2 


2 ^ 


TABLE V: The Gauge-Gauge-Scalar vertices. 


Vertices 

94>i4>2V 

Vertices 

50102 Vi 1/2 

Vertices 

50102^1^2 

H' 


hhW+W- 

4 

G^G'^WJZ^ 

ZEI 2 


^^2 —CK 

HHW+W- 

5^ 

4 

H^AW^Z, 

_(_,■ ffffz „2 

=n« 2 

G°G±IVT 

9 

2 

AAW+W- 

5^ 

4 

G^HWJZ, 

_99Z_a‘^ 

2 *W^/3-Q: 

H- 

2 

G°G°W+W- 

5^ 

4 

H^hW^Z^ 

_99^^2 

2 *VV"^/3-q; 

HH^W^ 

H- 

±*f S/3-a 

G+G-W+W- 

5^ 

2 

G^hW^Z^ 

_££Z 2 

2 *\4^*/3-a 

AH^W^ 

9 

2 

H+H-W+W- 

5^ 

2 

H^HW^Z, 

2 '^VF'^/3-Q: 

G+G-Z^ 

l^C2W 

hhZ^Zy 

5I 

8 

H^AW^A, 

=r^ 

^ 2 

H+H-Z^ 

i^C2W 

HHZ^Z, 

5I 

8 

G±G°W^Aj. 

“T 2 

hG^Z^ 

2 */3-a 

AAZy^Zi, 

5I 

8 

H^hW^A^ 

2 

hAZ^ 

2 ^/3-a 

G°G°Z^Z^ 

5I 

8 

G^HWifA^ 

2 '-'p —O' 

HG°Z^ 

2 ^/3-a 

G+G-Z^Z, 

al 2 

4 ‘'2W 

G+G-A^Z, 

e9ZC2W 

HAZ^ 

2 

H+H-Z^Z, 

al 2 

4 ‘'2W 

H+H-A^Z, 

e9zC2W 

G+G-A^ 

ie 

G+G-A^A, 


G^hW^A^ 

2 *P-a 

H+H-A^ 

ie 

H+H-A^A, 

e2 

H^HW^A, 

2 *^-0: 


TABLE VI: The Scalar-Scalar-Gauge and Scalar-Scalar-Gauge-Gauge type vertices and those coefficients 
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These coefficients are given by 


>'H+H-h = - [(2M^ - 2m'jj± - ml)sp-a + 2(M^ - ml) cot 2/3c/3-a] , 
Aaa/i = ^ [{2M^ - 2m\ - ml)sf}_a + 2(M2 - ml) cot 2/3c^_«] , 


^HHh — 


SjS-c 

2v L 


{2M^ — 2mjj — 'inl)sl_^ + 2(3M^ — 2mjj — ml) cot 2/3 s/3_q,C/3_c 


- (4M - 2mH - mh)cfi_, 


mt. 


m2 - 


mt 


'-Si3-aCj^_^ + 


m 2 - 


m 


h „3 


2v 


c^-Q,(cot/3 - tan/3), 


^hhh 2y ®/9—a T 

AcGh = - — Si3-a, 

^H±G^h = —ijnl - rn\±)ci3-a, 

V 

^AGh = --{jnl - m\)c^_a, 

V 

Xh+h~h ~ — 2(m2 — Tn\j) cot 2/3 s,3_q, + (2?tt.|^± + Tn\j — 2M2)c^_q 
^AAH — ~~ 2(M^ — mh) cot 2j3sf}—oi + {2m\ + mh — 2M^)c^_q, 


^HHH — — 


2v 

1 


2v 


2(m2 - mjj) cot 2/3s^_„ - 2{M^ - m]j)cp_c,sl_a + 


^GGH = — 

1 


m 


H 


2v 




^h±gth — -{mn - 'm^±)sp-a, 

V 

^AGH = -{m\j - m\)si3-a, 

V 


^Hhh = - 2yl\^2j3 “ 3M2)sin2a + M^ sin 2/3] , 


>^H±GTA = ±-("^A - fn]j±). 


The four point couplings are given by 


(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(AS) 

(A9) 

(AlO) 

(All) 

(A12) 

(A13) 

(A14) 

(A15) 

(A16) 

(A17) 


^H+H-AG = —-i^H+H-HS/3-a “ ^H+H-hC/3-a), (AlS) 

^G+G-AG = --(Ac+G-i^'S^-a “ ^G+G-hC/3-a) , (Al9) 

^AAAG = --{^AAHSjS-a “ ^AAhCjS-a), (A20) 

^AGGG = --(^GGHSp-a “ ^GGhCjS-a)- (A21) 
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Appendix B: Loop Functions 


The Passarino-Veltman functions 1071] are quite useful to systematically express the one-loop 
functions. First, we dehne A, B and C functions: 

d^k 1 


IOtt^ 




A-D _ 

(27r)^iVi’ 

4 - 1 ) / d^k [l,k^,k^k^ 


[C'o,C''",C^''](Pi,P2, (Pi +P2) ;mi,m2,m3) = p 


{ 2 tt)D N1N2 
4-D f [l,k^^,k^^k^ 


7 


(Bl) 

(B2) 

(B3) 


7(27r)^ N 1 N 2 N 3 ’ 

where D = 4 — 2e, and /r is a dimensionful parameter to keep the mass dimension four in the 
/c-integral. The propagators are dehned by 


Ni = k"^ - ml + ie, N 2 = {k + pi)^ -ml + ie, N 3 = {k + pi + p 2 f‘ -ml + is. 


(B4) 


The vector and the tensor functions for B and C are expressed in terms of the following scalar 
functions: 


(B5) 

(B6) 

C^=p^Cii+p^Ci2, (B7) 

c^'' = p‘tp'iC 2 i + P^ 2 P 2 C 22 + Kp2 " + PiP‘i)C 2 ^ + 9^''C2,. (B8) 


By counting the mass demension of the above functions, we can find that the divergent part is 
contained in A, Bq, Bi, B 21 , B 22 and C 2 a- All the scalar functions are expressed by the divergent 
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part and finite part as 


A{m) = (A + 1 — Inm^) , 

Bo = A — f (IxIuAb, 

Jo 

A 

Bi = -—+ dx{l - x)lnAB, 
^ Jo 

/ dx{l — x)^ In As, 
Jo 


B21 = -- 
3 

n ^ f 2 , 2 /\ f 2 , 2 

^22 = ^ I mi + m2 - Y 1 A + - I mi + m2 - Y 


dxAB In As, 


^0 = - 
Cii = - 

Cl 2 = — 
C 2 I = — 
C 22 = — 
C 23 = — 


[ dx [ 

Jo Jo 

[ dx [ 

Jo Jo 


dy 


y 


Ac’ 

y{xy - 1) 


dy- . 

0 Ac 


/■‘dx 

Jo Jo 


[ dx [ 

Jo Jo 

[ dx [ 

Jo Jo 

[ dx [ 

Jo Jo 


dy- . 

0 Ac 

y{l - xyY 


dy- 


dy 


dy 


Ac 

y(i - yf 


A 


c 


y{l - xy){l - y) 


C 2 A = ^ -I 


[ dx [ 

Jo Jo 


Ac 

dyy In Ac, 


where 


Ab = —x{l — x)p^ + xm\ + (1 — x)m 2 , 

Ac = y'^iPix + P 2 f + y[x{pl -q^ + ml- ml) + m^ - ml - pi] + ml, 


(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(BIT) 

(B18) 

(B19) 

(B20) 

(B21) 

(B22) 


and the divergent part A is given by 


A =- 7 c + In dvr + In 


with 7 g being the Euler constant. It is convenient to define the following functions 


(B23) 


10811 : 


B2(/,mi,m2) = B2i(/,mi,m2), (B24) 

B3(p^,mi,m2) = -Bi{p^,mi,m 2 ) - B 2 i(p^, mi, m 2 ), (B25) 

B4(p^,mi,m2) = -mlBi{p^,m2,mi) - mlBi{p^,mi,m2), (B26) 

B5(p^, mi, m 2 ) = x4(mi) + A(m 2 ) - 4^22(p^, mi, m 2 ). (B27) 
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Appendix C: IPI diagrams 


In this section, we give the analytic expressions for the IPI diagram contributions to one, 
two and three point functions by using the Passarino-Veltman functions defined in the previous 
section. We calculate IPI diagrams in the t’ Hooft-Feynman gauge in which the masses of Nambu- 
Goldstone bosons mQ± and niQO and those of Fadeev-Popov ghosts and rricz same as 

corresponding masses of the gauge bosons; i.e., mQ± = = my/ and m^o = ruc^ = mz- IPI 

diagrams with bosonic external lines are separately calculated by the fermion-loop and boson-loop 
contritbutions. We denote the fermionic- and bosonic-loop contributions by the subscript of F and 
B, respec tively . Throughout this section, we use the shortened notation of the Passarino-Veltman 
functions 1071] as 


A(V) = 


IGvr^ 


A{mx) 




{p^;X,Y) = 


16vr2 

2 . 


Ci,ij{X,Y,Z) = Y^Ci^ijiPi,P2,{Pi +P 2 ) ■,'rnx,mY,mz). 


(Cl) 

(C2) 

(C3) 


1. One-point functions 


The IPI tadpole diagrams for h and H are calculated by 

= - E 

f 


(C4) 

(C5) 


T^lPI _ „ 
^h,B — 'S/3-C 


3gmwA{W) -h -gzmzA{Z) - 2gm^ - gzm^z 


— - ^AAhA{A) - \HHhA{H) - 3XhhhA{h) 

— ^G+G-hMG"^) ~ ■^G0G0/i^(G'°), 


tIPI - r 

2 H,B — Cp-a 


3gmwA{W) + -gzmzA{Z) - 2gmw - gzm% 


- ^h+h-h-^{B^) - ^aahA{A) - 3\hhhA{H) - XHhhA{h) 

— ^G+G-HAiC^) — Xq0q0hA{G^). 


(C6) 


(C7) 
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2. Two-point functions 


The IPI diagram contributions to the scalar boson two point functions are calculated as 


f 

/ 

f 


Aif)+l2mj-^ Boip^-JJ) 


.2 P 




Aif) + [2mj-i-]Bo{p^-J,f) 


f 


Aif)-^Boip^-JJ) 


nJS(pV = -E 


AmlNi 




A{f)-^Bo{p^;f,f) 


(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 


'^}KHp^)b = sin2(/3 - a)(3m^ - p^)Bo{p^] PP, IT) + y [4 - siii^ijS - a)] A(W) 


+ Y - P^)Bo{p^] Z, Z) + ^ [4 - sin2(/3 - a)] A{Z) 

^2 


al 


- cos 2(/3 - a) [2^(PT) - A{H^) + (2m|^± - + 2p^)Bo{p^-W, H^)] 


- y cos^(/3 - a) [2^(Z) - A{A) + (2m^ - mz + ‘^P^)Bo{p‘^; Z,A)] 


sin^(/3 — a) + 


1 


{2g^m^ + g\m\), 

- 2Xfj+fj-f^fiA{H^) - 2\AAhhA{A) - 2\HHhhA{H) - l2\hhhhA{h) 

— 2\q+q-i^j^A{G^) — 2\qQqQi^j^A{G^) 


+ Xl^H-hMp"-, + XI+c-hMp'-, G±) + 2A|,+G-h^o(p^ 

+ a, A) + 2XQOQOfiBo{p^-, G°, G°) + X^QO^Baip"^A, G°) 


+ ‘^AHHh^oiP^'i H, H) + 18A|;j;ji?o(p^; h, h) + 4X^fjf^f^Bo{p'^; h, H), (C13) 
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cos^(/3 - a){3m‘^ - p^)Bo{p^-, W,W) + ^[4.- cos^(/3 - a)\A{W) 

+ - ot){3m\ -p^)Bo{p^;Z,Z) + ^ [4- cos^(/3 - a)] A{Z) 

q 2 

- ^ sm\f3 - a) [2A{W) - A{H^) + (2m^± - + 2 p^)Bo{p^; W, H^)] 

- ^ sin2(/3 - a) [2A{Z) - A{A) + {2m\ - m| + 2p^)B^{p^-Z, A)] 




{ 2 g^m^ + glm%), 


- 2Xh+jj-jjjjA{H^) - 2XaahhA{A) - 12XhhhhA{H) - 2XHHhhA{h) 


- 2Xq+q-jjhA{G^) — 2Xqoqohh^{G'^) 


+ xjj+H^HBo{p^-, + a|+g-h^o(p"; G±) + 2A^+g-h^o(p"; H^,G^) 

+ ‘^^AAH^oiP^'tA, A) + 2X^QaQ0jjBQ{p^] G°, G°) + X\qojjBq{p^] A, G°) 

+ H, H) + 2A^;j^i?o(p^; h, h) + (C14) 


n/Fh(p^)B — Sp-aCp_a 

X {<7'(3m2^ _ py) _ 

+ “ P^)Bo{p^] z, Z) - ^A{Z) 

+ ^[ 2 A{W) - A{H^) + (2m^± - + 2 p^)Bo{p^-W, H^)j 

+ ^[274(Z) — ^(A) + (2mA ~ ''^1 + ‘^P^)Bq(p^] Z, A)] — {2g‘^m^ + g\m%)^ 

— ^H+H-Hh^{B^) - XAAHhM^) - ^^HHHhA(H) - 3XHhhhA{h) 

— ^G+G-Hh^(G^) — ^G°G°Hh^{G^) 

+ ^H+H-h^H+H-nBoiP^', H^, H^) + ^G+G-h^G+G- hBo(p^', G^, G^) 

+ ‘^^H+G-h^H+G-HBo{p‘^', H^,G^) 

+ ‘^^AAh^AAHBo{p‘^] A, A) + 2XfiQ0Q0 XqocohB o{p‘^G^, G^) 

+ ^AG°h^AG°HBo(p^]A, G°) + GXnHh^HHHBoip^] H, H) 

+ Q^hhh^HhhBoip'^; h, h) + 4XHhh^HHhBo{p‘^; H, h), (C15) 
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^^aa{p^)b — + 9z^{Z) — -{2g^m\i + gzrn%) 

- ^ [2AiW) - A{H^) + (2m^± - + 2p^)Bo{p^-, W, H^)] 

- ^ cos^(/3 - a) [2A{Z) - A{h) + (2m| - m| + 2p'^)Bo{p^; Z, h)] 

- ^ sin2(^ - a) [2A(Z) - A{H) + {2m\ - m| + 2p^)Bq{p^-Z, H)] , 

- 2Xh+h-aaA{H^) - 12XaaaaA{A) - 2XaahhA{H) - 2XAAhhAW 

- 2Xq+q-aaA{G^) — 2Xaag^g^A{G^) 

+ 2\XH+G-A?Bo{p^;H^,G^) + AX\AhBo{p'\A,h) 

+ ^^\ahBo{j)^', A, H) + X\qqj^Bq(j?] h, G^) + X\qq^Bq(ji^-, H, G^), (C16) 


B^Ag{P^)b — S/S-aCjS-a 

X [2A{Z) - A{H) + {2ml " "^1 + V)i^o(p'; H)] 

- ^[^A{Z) - A{h) + {2ml “ "^1 + 2p^)^o(p^; Z, /i)]| 

“ ^H+H-AG°A{H^) — 3XaaAGoA{A) — XagohhA{H) — XAG°hhA{h) 

- ^g+g-agoA{G^) - 3XagogogoA{G^) 

+ ‘^^AAh^AG°hBo{p^'i A, h) + 2Xaah^ag°hBo{p^', A, H) 

+ ‘^^AG°h^G°G°hBo{p‘^] G^, h) + 2Xag°H^G°G°hBo{p^] i H)- (C17) 


The Z-A mixing is given by 


2 mj 


Aza{p^)f = ^ -^'rnzNi^fBo{p‘^; /, /), 
/ ^ 

B.za{pIb = m-z 


(C18) 


+ Bo){p^;A, H) - + Bo){p^;A, h) 


- ^^c^_„(2i?i + Bq){p^-G l H) - ^A^sp.^{2Bi + Bo){p^-, G^ h) 

2 2 

- ^-Sfj-aCp_a{Bi - Bq){p^- H, Z) + ^Sf}-aCp_a{Bi - Bo){p^- h, Z) 


(C19) 


The G.I. part appearing in Eq. (jllSp is given by 


^za{pHqj = Aza{p‘^)f 
2mz 


+ 


^AAFiSfj-a{‘2'Bi + Bo){p ; A, H) — XAAhCf}-ai‘^Bi + Bq){p ',A,h) . (C20) 
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The IPI diagram contributions to the gauge boson two point functions are calculated as 


KwiP^)F = (2p^Bs - 54) (p"; /, /'), 

fj' 

= E 8 e^Q}Nip^Bs{p^- f, /), 

/ 

= E \2p\2ifQf - islrQ})Bs\ (p 2 ; /, /), 


(C21) 

(C22) 

(C23) 


nFi(p")F = E aW [ 2 p\As^^Q) - As^wQflf + - 2 l}fBo\ (p^ /, /), (C24) 

/ 


nw?vy(p^)s = 9^|^-B5 (p^;^,-H'^) + ^sin2(/3 - a)B^{p^\H,H^) 
+ I cos^(/3 - a)B^{jP‘; h, H^) 


+ sin (/3 - a) + -B 5 j (p ; h, W) 

+ cos^(/3 - a) (^m^Bo + ^^5^ (p^; F, IT) 


+ 


^ ) B^ + {rriy^ — 4s^ym^^ + — 8p c^)i?o 


(p2;Z,lT) 


+ 2s^ i?5 + (2m^ — 4p^)i?o (p^; 0, W) — —p^ 


= e^B5{p^-,H^,H^) - eV 


12^3 +5^0 (p^bh,Vh) + 


B%\p‘')b = ^B,{p^-, H^, H^) - egzp^ ( IOB 3 + + ^ ) (p^; W, IT) 




civ 

nFJ(p")B = ffl 11^5 (p'; + l sin2(/3 - a)i?5(p'; ^) 

+ ^ cos^(/3 - a)B^{p^] h, A)] 

+ sin2(/3 - a) ^m|So + ^-85^ (p^; h, Z) 

+ cos^(/3 - a) (mlBo + ^B^'] (p^; H, Z) 


+ 


OQ 1 O 

(2m^ - -jP‘^)Bo - Op^Ss] (p^; IT, IT) - -p^ 


(C25) 

(C26) 

(C27) 


2 sw 


W,\p^)b - 

Ct, 


’-VI/ 

where the fermion-loop contributions are the same as those in the SM. 


(C28) 
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The fermion two point functions can be decomposed into the following three parts 




tIPI 


IPI 


tIPI 


(C29) 


Each part is caluclated as 


~ + l)(p^; /) 7) “ 9z{v^f + Oj)( 2 i?i + l)(p^; /, Z) 

■,2 


-^(2i?i + l)(p2;/',VE) 


m. 


idfB^ip^- f, h) + {euYBiip^- f, H) + ijBiip^-, f, A) + Bi{p^-, f, GY 

2e2 I ^2 t2 


:f ^2 i 


rnl + rnl, . 9 , + ^ , ,, 


= -2<7|p/a/(2Si + l)(p^ /, Z) - ^-{2B, + l)(p^ /', IT) 


■,2 


2. 


+ 


rrYf — m% „ , , w-fC? — 

f ^- 5 i(p 2 ; /', G±) + f f f f BYp^-, r, H^), 


n//,s(p') = -2e^Qj{2Bo - l)(p";/,7) - ^gUvj - «/)(2i?o - l)(p";/, Z) 


. 2^2 


,2 /„,2 


m 


+ 


i^YBoipY /, /i) + {enYB^ip^- f, H) - YjBoip^-, f, A) - /, 


^2 


2 . ^ 


- 2 - 


m 


f 


[So(p'; /', G±) + Y^pBoip^-, /', , 


(C30) 


where u/ and af are the coefficient of the vector coupling and axial vector coupling of Zff vertex 
given as 


~ 2 ^wQf^ ®/ ~ 2 ■ 


I 


f 


(C31) 


3. Three-point functions 


In this subsection, we give analytic expressions for the IPI diagram contributions to the three 
point functions. The assignment for external momentum is taken in such a way that pi and (^ 2 ) is 
the incoming momnetum of h (h), V (V) and / (/) for the hhh, hVV and hff vertices, respectively, 
and q = Pi + P 2 is the outgoing momentum of h for all the above vertices. 

First, the IPI diagrams for the hhh coupling is calculated as 


_ /V'i* 

^hhh\Pl^P2^Q )F — — 2_^ 3 Is/j; 

/ 


Bo{pY f, f) + Bo{pI, f, f) + Bo{q^, /, /) 


+ (4m^ -q'^+Pi- P2)Co{f, /, /) 


(C32) 
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= ^rn^sp-a lQCo{W,W,W) - Cq{c^, c^) 


,3 „3 


± „± „±^ 


- -l^BawSfi-a 

+ f ™l4-« 

^2 




SVV tTj± 


16Co(^, Z, Z) — Cq{cz.,cz.,cz) 


9z^z 


Sj3-c 


+ ^^G+G-hS0-aGhhh(^^G'^,G^) + ^^H+H-hC0-aGi 


hhh 


Z, zyy^cl^a^ {A, Z, Z) 

{W,H^,H^) 


+ yAH+G-/^5/3-aC^-a[QV/^^(V^^, G^, H^) + H^, G±)] 

+ ^^G°G°hSp-aGhhh (-^J G^, G°) + ^^AAhC^p-aGhhh {^■> ^) 


9% 


+ ^AyiG'0/iS/3-aC^-a[<^hWi* G'‘^) + G^^h (^) ^)] 

+ 2ff^”^VKS/j-a[So(p?, + Bo{pI, W, W) + Bo{q^, W, W)] - 3g^mws^-a 
+ 9z^zSfj-a[Bo{pl, Z, Z) + Bq{p% Z, Z) + 5o(g^ Z, Z)] - ^9%'mzsp_a 
+ 2\H+H-h^H+H-hh[BM, + Bo{pI + Bo{q\ H^,H^)] 

+ “^^hG+G-^hhG+G- [Bo{Pi, G^) + -Bo(P 2 ) G^, G^) + Bo{q‘^, G^, G^)] 

+ ^^H+G-h^H+G-hh [Bo {pI ,H^,G^) + Bo{pIH^,G^)+ Bo (q^ ,H^,G^)] 


+ ‘i:2iAAh2^AAhh[Bo{p‘i, A, A) + Bo{p 2 , A, A) + Bo{q^, A, A)] 

+ '^^G°G°h^G°G°hh[Bo{p\-, G°, G°) + i?o(P 2 ) G°, G°) + Bo{q^, G°, G'^)] 
+ ‘2^AG°h^AG°hh[Bo{Pi, G°) + Bo{p 2 , A, G°) + Bo{q^, A, G°)] 


+ ^^HHh\HHhh[Bo{pl H, H) + Bo{pI H, H) + Bo{q\ H, H)] 

+ ^2XHhh^Hhhh[Bo{pi, h, H) + i?o(P 2 ) h, H) + Bo{q^, h, H)] 

+ '^ 2 \hhhKhhh[Bo{pl, h, h) + Bo{pI, h, h) + Bo{q^, h, h)] 

- 2Xj,+j,.^Go{H^,H^,H^) - 2A|+c_^Go(G±,G±,G±) - 8 A^ogo/,Go(G°, G^, G^) 
“ ^^AAhGo{X^, X\.) — 8X%fjf^Go{H, H, H) — 2\QX\y^^Go{h, h, h) 
-2AH+H-/^A^+G-h[^o(G±,F±,//±) + Go(//±,G±,F±) + Go(F±,//±,G±)] 
-2AG+G-hAi+G-h[C'o(^^,G±,G±) + Go(G±,F±,Ty) + Go(G±,G±,//±)] 

~ 2 XAAh 2 ^AG°h[Go{G^, A, A) + Co{A, G*^, A) + Go{A, A, G^)] 

- 2AGOGO/.A^GOh[C'o(Al, G°, G°) + Go(G°, Al, G°) + Go(G°, G°, Al)] 

- 8 XHHhXHhh[Go{h, H, H) + Co{H, H, h) + Go(F, h, H)] 


- 24XhhhXjjhh[Go{h, h, H) + Go{H, h, h) + Co{K H, h)], 


(C33) 
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where 


^hhl^ (^! — P 1 C 2 I + P 2 C 22 + 2piP2C'23 + 46*24 ~ 7 : ~ (^ + Pl)(PlC'll + P 2 C 12 ) + Q'PlC'o (^j Y, Z) 


+ 


+ 


p\C 2 l + pIC 22 + 2 piP 2 C '23 + 46*24 “ 77 + ( 3 pi “ P2)(Pl6ii + P2C12) + 2 pi{pi - ^2)69 {Z, X, Y) 


P1C21 + P2C22 + 2P1P2C23 + 4624 “ 77 + ( 3 pi + 4p2)(piC*ii + P2C12) + 2 q{q + ^2)6*0 {Y, Z, X), 


Cj:,f{X,Y,Z)^ 


1 


PiC*2i + P2C22 + 2P1P2C23 + 46*24 “ 2 ^^'2)(pi6ii +P2612) + 4 pi • qCQ (X, Y, Z) 

+ P1C21 + P2C22 + 2P1P2C23 + 4624 “ 2 ^^'2(^1611 +P2612) — pi{pi + 2 p 2 )Co {Z, X,Y) 

+ Pi6*21 + P2C22 + 2P1P2C23 + 4624 ~ 2~ 2^2(^1611 + P2C12) — q{pi — P2 )Cq {Y, Z, X). 

The hff vertex can be decomposed into the following 8 form factors 

r}^Mplplq^) = 

^hff + + 1^2 T)l// + i>l + 1^2 I^Fhff + 1^1 ^2 F^ff + ]/>l i>2 l^FhJf ■ 

Each form factor can be calculated by 

Fhff = “ 2 ( 7 |u^(uj — aj)sp-aCo{Z, /, Z) 

~ + Pi (6*11 + 621) + P2 (C'i 2 + 6*22) +Pi • ^2(2623 — Co ) + 4624 — !](/) 7 ) /) 

+ dzi'^'f ~ Q/)[^/C*o + Pi(6ii + 6*21) + ^2(^12 + 6*22) + Pi • P 2 ( 26*23 ~ Co ) + 4624 — !](/, Z . 

{d)^C^ff{f, K f) + {Ci)^C^ff{f, H, f) - C^ffif, 6°, /) - ffCfffif, /) 


(C 34 ) 


(C 35 ) 


n y2 


-c, 


^7T? ^ r 


h y2 


CCfif, C^,f) + n 


FSFffi TT± 


nil 

V 


f{ml?KhhCo{h, /, h) + 2{^{,)^XHHhCo{H, /, H) + 2CiciXHhh[Co{h, /, H) + Co{H, /, h)] 
- 2Agogo/.6o(GO, /, G°) - 2 efXAAhCoiA, f, - C/AuGOh[Go(7l, /, G°) + Go(G°, /, ^)]} 
+ AG+G-/^Go(G±, /', G±) + 

+ l>^H+G-h{Cf + C/')[Go(G±, /', + CoiH^J', G±)]} 


-S/3-C 


CXffiW, f, G±) + C^J/ (G±, /', fE) 


hff 


- ^-^ifCp-o. [clffiW, /', i6±) + C^^f/{H^,f\ fE) 




S/3-C 


/, + C^f/(C^, f, Z) 


SFViriO 


9z, 


f,A)+ C^ff' {A, /, Z) 


SFV/ 


(C 36 ) 
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.2 


-1 


Fhff = -Cf)[CoiG^,f,H^) 


-Co{H^,f,G^)] 


9 

T"^- 

9 


'clff{W, /', G±) - /', VT) 

^^fcp-a [c}^f7iW, f, H^) - /', 

9lvfIfSf,_^ [GXf/iZ, /, GO) - Gf//(GO, /, Z)' 

[Glff{Z, f, - G^f/iA, /, Z)' 


(C37) 


j,vi _ 2"^L/ 

9 f' f 

A _ LeJ I 


gUvf + a/)(Go + 2Gii)(/, Z, /) + (Go + 2Gii)(/, 7 , /) 


+ 5'^a(Go + 2Gn)(/',M^,/') 


- sp_^g%v{v} + a}){Go + Gii)(Z, /, Z) - + Gii){W, /', 

(e()2(Go + 2Gii)(/, /i, /) + (Ci,)'(Go + 2Gn)(/, F, /) 

+ (Go + 2Gii)(/, GO, /) + e^(Go + 2Gn)(/, 7l, /) 

(m} + m2,)(Go + 2Gii)(/', G±, /') + (m^ef + m},er)iCo + 2Gii)(/', /') 

{6(C()"A;,;,;,(Go + Gn)(/i, /, /i) + 2(e^)2AH/fh(Go + Gn)(77, /, 77) 

+ 2^l^{,XHhh[{Co + Gn){H, /, h) + (Go + Gii)(/i, /, if)] 

+ 2\GOGOh{Co + Gii)(GO, /, GO) + 2e|A^Ah(Go + Gii)(Al, /, Al) 

+ C/Aago/^[(Go + Gn)(A, /, GO) + (Go + Gii)(GO, /,. 


4 

rUl} 


in r/ r/ 

yi 


m 


f 


A, 


G+G-/1 / 2 


(mj + mj,){Go + Gn)(G±,G±) - 


^H+H- 


^H^G-h i_2t , _2 t _\ur^ , ^ ^/^± ft tj± 


,± ,, ^"■^■'^(m}C^ + m2,^2^)(Go + Gn)(ii±,/',ii±) 


m). 


^2 (^K/ + ^K/')[(^o + Cii)(G±,/',ii±) + (Go + Gii)(ii±,/',G±)] 

2 

\fi-a{2Go + Gii){W, /', G±) + S;3-a(-Go + Gii)(G±, /', W") 

- e/'C^-a(2Go + Gn)(ty, /', ii±) - irCp_^{-Go + Gn)(ii±, /', H^) 


4u 


2 "^/ 


S;3-«(2Go + Gn)(Z, /, GO) + S;3-a(-Go + Gn)(GO, /, Z) 

- ^fC0_M + ^^ii)(^, /, Al) - ifCp_^{-Go + Gii)(AI, /, Z) 


(C38) 
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f,V 2 _ 2 "^/^/ 
^hff - 


V 


-fj 


gl{v) + aj){Co + 2Ci2)(/, Z, /) + e^Qj{Co + 2 Ci 2 )(/, 7 , /) 


2/^2 , 


+ 9^'^Cl!iCo + 2Cu){f,Wj') 


- S 0 -ag^viv] + a})Ci 2 {Z, f, Z) - sp_^^vCi 2 {W, /', VT) 




h y3 


ft ml, 

+ 4^ 


{(i)HCo + 2Ci2)(/, h, f) + (^ifiCo + 2Ci2)(/, H, f) 

+ (Co + 2 Ci 2)(/, G°, /) + C^(Co + 2 Ci2)(/, f) 

(mj + nij,){Co + 2Cl2)(/^ G^, /^) + (mj^'j + + 2Gl2)(/^ f') 


-:J-{mif^hhhCi2{h, f, h) + 2{i^^fXHHhCi2{H, f, H) + 2e^C:^AH/.h[Gi2(i/, /, /i) + Gi2(/i, /, //)] 

+ 2Agogo„Gi2(G°, /, G°) + 2ef\AAhCi2{A, /, A) + 2e/A^GOh[f^i2(G°, /, A) + Gi2(A, /, G°)]} 
(my + my,)Gi2(G=^, /', G=^) - + "i/'^/')C'i2(^^'^, f, H^) 


^G+G-ht 2 


+ m2,eyO[Gi2(G±, /', if±) + G±)] 

s^_„(2Go + Gi2)(Vr, /', G±) + s^_„(-Go + Gi2)(G±, /', P^) 


V 

4 ?; - 


2 m^.. r 


4!!!/ 

8 V 


- e/'C^_„(2Go + Gi 2)(PP^, /', - C/m^_„(-Go + Gi2)(i?±, /', PF) 

s^_„(2Go + Gi2)(Z, /, GO) + s^_„(Gi2 - Go)(GO, /, Z) 

+ ^/C/3-a(2Go + Gi2)(Z, /, A) + ‘^/C^_„(Gi2 — Co){A, f, Z) 


(C39) 


TTh^ TTh^ 

F^ff = + 2Gn)(/, Z, /) - g^^C({Co + 2Gn)(/', PP^, /') 


2u 


+ 2s^_„4t-yayi;(Go + Gii)(Z, /, Z) + s^_^^z;(Go + Gii)(Py, /', PP^) 


m 


+ ..3 




(m} - mj,)iCo + 2Gii)(/', G±, /') + (m^^/ - m^,4)(Go + 2Gii)(/', /') 


f/ [Tit 


'G+G-/1 


{mj - mj,){Co + Gii)(G±, /',G±) - 


^H+H-h / 2 c2 2 c2 


{mjef - m},4)(Go + Gn)(i7±, /', 


yZ J J yZ 


- m}tCf,)[{Co + Gn)(G±,/',/7±) + (Go + Gii)(F±,/', G±)] 
S;3-a(2Go + Gn)(PP^, /', G±) + s^_„(-Go + Gn)(G±, /', PT) 


+ 


V 

4 V 


- j2Go + Gn)(PP^, /', - ^f'C^-ai-Co + C'ii)(i7±, /', PP^) 


TTlj 


L 


+ gpfVf^ s/3-„(2Go + Gn)(Z,/,G^) + s/3_„(-Go + Gn)(GO,/,Z) 


+ CfCis-a{2Co + Gii)(Z, /, ^) + ifCp_^{—Co + Cii){A, f, Z) 


(C40) 
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TTL^ Tirt^ 

^hff = -i^l9lvfaf^iCo + 2Ci2)if,Z,f)-Cl'g^^iCo + 2Ci2)if,W,f) 


2v 


+ 2si3.^g^zyfafvCi2{Z, f, Z) + s/s-a^vC^iW, /', VT) 


{m} - m}){Co + 2 Ci 2 )(/',G±,/') + (m}?/ - m),i}){Co + 2Cu)if, , f) 

{mj - m‘j,)Ci 2 {G^, /', G^) - f, H^) 


+ 


- m},Cf,)[Gi 2 iG^,f, H^) + Gi 2 {H^,f, G±)] 

m), r 


4 V - 


s^_„(2C7o + Gi2){W, /', G±) + Sfi-o^i-Go + Ci2)(G±, /', 

- irCp_^{2Go + Gi2){W, /', H^) - ^f,c^_^{-Go + G^m^J', W) 


rrij 


+ gPfVf^ [s^-a(2Co + Gi2 ){Z, f, G^) + S0.a{-Co + Ci2)(G^ /, Z) 
+ ?/C/3-a(2C'o + Gi2 ){Z, /, + ^fCj^_a( — Go + C'l2)(^, /, 


(C41) 


o 

fTTif r 

Z?-^ _ CJ _L 

^ y ^hff - ^2 




(C{)2(Cn - Gu){f,hJ) + (Ci,)'(Cii - Gi2)if,H,f) 


- {Gn - Gu)if,G^,f)-ej{Gn - Gi2){f,AJ) 

\cu-Gi 2 ){f',G^,f)+^f^f,iGn-Gu)if,H^,f) 


_G'GGl 

^2 


r 

4 


5| ^ 


S0.ai-2Go - 2Gn + C'i 2 )(VF, /', G±) + s^.^{-Go - Cn + 2 C'i 2 )(G±, /', 14^) 

+ e/c^-a(-2C'o - 2Cn + C'i2)(Vr, /', if±) + CfC^_^{-Go - Gu + 2 Gi 2 ){H^, /', 
s^_^{- 2 Go - 2Cn + Ci2)(Z, /, G^) + s^_„(-Co - Cn + 2 Ci2)(G0, /, Z) 


+ C/c^-a(~2Go — 2Gii + C'i2)(^, f,A) + ^/c^_Q,(—Go — Gn + 2Gi2)(^, /, 


(C42) 


^ PT _ 9 ^ 


(?) 


s^_„(2Go + 2Gii - Gi 2)(W^, /', G±) - s^_„(Go + Gn - 2Gu){G^, /', 

- e/c^_„(-2Go - 2Gn + Gi2)(44^, /', JGq + Gn - 2 Gi 2 )(F±, /', W^)' 

- gllfVf L_„(-2Go - 2Gn + Gi 2 )(^, /, G^) + s^_„(Go + Gn - 2 Gi 2 )(G? /, Z) 


+ 2Go — 2Gn + Gi 2 )(Z, f,A)+ C/C/ 3 _a(Go + Gn — 2Gi2)(^, /, Z) 


(C43) 
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where 




[mf'C'o + PiiCii + C 21 ) + ^ 2(^12 + C 22 ) + 2pi • p2iCi2 + ^^ 23 ) + 4C'24](X, Y, Z) — 


Clff{X,Y,Z)^ 


'hff 


[p\{ 2 Cq + SCii + C21) + ^2(2^12 + 6*22) + 2 pi ■ p 2 { 2 Cq + 2 Cii + C12 + (^23) + 4 ( 724 ](X ,~ 


{X, Y, Z) = \pi{C 21 - Co) + p^(C 22 - C 12 ) + 2 pi • P2{C23 - Cu) + 4C24](X, y, Z) - 


'hff 


(C44) 


The IPI diagram contributions to the form factors of the hZZ and hWW vertices which are 
defined in Eq. (I125P are calculated as 


T^ 1 , 1 PI/ 2 2 2 \ 

^hzz \Pi^P2^ Q )f — / , 

/ 


lQm?jm\Ni 


{vj + a}) Boip'i, f, /) + Bo{p^, f, f) + 2 BQ{q\ f, f) 


+ ( 4 m 2 -p\- p^Coif, /, /) - 8 C 24 (/, /, /) 

- (vj - aj) Bo{pI, /, /) + Bo{pI, /, /) + (4m^ - q^)Co{f, f, f) 


t^ 2 ,ipi^ 2 2 2 ^ 32m^m|A^/ 

^hzziPi^P2^^ )f = - 2 ^- ^^^3 - 

/ 

i^f + fl/)(4(723 + 3(7 i2 + Cii + Cq) + (vj — ajr){Ci2 — Cu) (/,/,/), 

t.3,1PI, 2 2 2 ^ ^64m2m|A^/ , r. 

^hZZ iPi^P2^Q )f = 2_^ -^-3-+ C12 + Co){f, f, /), 


(C45) 


(C46) 

(C47) 


T^hiPi . 2 ^2 2 ^ _ Y- at/ 


/./' 


^Boipl /, /') + Bo{q\ /, /) + -Boipl /, /') 


- 4 C 24 (pi,pi, 97 f, f, f) + + 2m^, - PI - pl)Co{f, /', /) 


+ (m/ o- ruf/), 


(C48) 


—Arriu/rriiN/ 


^IwUpIpI 9')f =- ^3 ^ ' (4C23 + 3Ci 2 + Cii + Co) (/, /', /) + (m^ o myO, (C49) 


—AmfirmlNi 

{pi,pi, 9^)f =- 3 -(Cii + C 12 + Co) (/, /', /) + (m/ o m//), 


■^ 3 , 1 PI ^ 2 2 2 ^ 

/iVPVP 


(C50) 
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= 2gl\G+G-^mlrS^wCo{G^,W,G^) 

+ 9^ fBW S p-a\^^G]!^yyi{W, W, W) — 2<^C2 a{c^,C^, C^) + S^C^yyi{G^, W, W) + S^G'^yyi{W, W, G^) 
_2£fm^s^_„Go(Ty,G±,H^) - (^ - 4)f|l[G24(W, G±, G±) + C24(G±, G±, W)]} 

Cy/ Cy, J 

3 

H—— 2m^ ^s‘jj_^Go{Z, h, Z) + c^_^Gq{Z, H, Z)'j + s^_„[C'24(G^, h, Z) + G24(^, h, G^)] 

+ 4-a [C2 a{A, h, Z) + G24(^, h, A) + C2a{G\H, Z) + G24(^, H, GO) - G24(A, H, Z) - G2a{Z, H, A)] } 

+ 2g|m^|3A/i/i/iS^_QGo(/i, Z, h) + XhH hC^-aGoiH, Z, H) + A/f/i/iS^_„c^_„[Go(-ff, Z, h) + Go{h, Z, 


- 2g|(4 - 4)2 [AG+G-h^^24(G±,G±,G±) + A^+^-;,G24(i^^, i^^)] 

“ ‘^9z^‘p-a 3A/j/i/jG24(^, G^, h) + XHHhC24:{H-, A, H) + Agg/iC'24(G0, h, G^) + A2l2l/iC'24(^) -ff) ^) 
“ ‘^9z^‘0-a ^^hhhC24{h, A, h) + XHHhC24:{H, G^, H) + XAAhC2A{A, h, A) + XGGhC24{G^, H, G®) 

- 2g|s^_„c^_„AH/^h[C'24(/^, GO, i?) + C24{H, G^ h) - G2A{h, Al, H) - G24{H, Al, h)] 

- 2glsp_^Cp_^XAGh[G24{A, h, G^) + G24{G^, h, A) - G24{A, H, G^) - G24(G0, H, Al)] 

2 2 

+ YAG+G-h(4 - 4)'^o(<?^G±,G±) + ^Ah+h-/x(4 - 

+ ^^GGhBoiq^, GO, GO) + ^XAAhBoiq'^, A, A) + ^XuHhBoiq^, H, H) + -^XhhhBo{q^, h, h) 


A 

,3 *VK„ 


- g'^^mws^_a[Bo{pi,W,G^) + Bo{pi,G^,W)] - ^mzS0_jBo{p(,h, Z) + Bo{p^,h, Z)] 






-w 


- Qg^CyrmwSi3_aBo{q'^, W, W) + 4.g^c^mwSi3_a, 


(C51) 
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idziT^z) ^'^hZ^ziPl^Pl^ = ‘^g'mwCwSp-aChVV2{W^ W) - 2gcwmwSf}-aCi223{c^, c^,c^) 
+ d'lBWSw<^W^l3-a[ChVV2iG^^ + ChVV2i^^ 

- gmwiclr - sl,)sUC^^V2iG^,G^,W) + C}^^^2iW, G±, G±)] 

+ ^iBz[Ghvv2i^^^^G^) + C](yy 2{G^,h, Z)] 

+ ^'>BzS^-a[GhVV2i^^ h,G^) + G§yV2{G^j h, Z)] 

+ ^1BzS(3-aC^-a[GhVV2{^j ^) + GXvV2{^^ H, G°) — C^yV2{^^ G[, A) 

+ G^yV 2{^: h, Z) + G^yV2{G^, H, Z) — C^YV2{^-I ^1 ^)] 


- 2(C^ - 4)2 [AG+G-hCl223(G±, G±, G±) + Ah+//-/^Ci223(^^, 


- 2s 


2 

p—a 


3^hhhGi223{h, G^, h) + XHHhGi223{H, -f^) + AGGhC'l223(G°, h, G^) + AaA/iC'i223(^) 


- 2c' 


73—a 


^XhhhGi223{h, A, h) + XHHhGi223{H, G°, -ff) + A^AhC'l223(^) + AGGhC'i223(G°, G°) 


- 2s^_qC^_q,Ah/i/i[C'i223(/i, G^,H) + Gi223{H, G°, /l) - Gl223{h, A, H) - Gi 223{H, A, h)] 

- 2Sf,_^Cp_^XAGh[Gl223{A, h, G^) + Ci223(G0, h, Al) - Ci223(^, H, G^) - Ci223(G°, H, Al)], 


p 3 ,lPI 
^ hZZ 


ipl,pl,q‘^)B = 0, 


(C52) 

(C53) 
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T^l.lPI ! 2 2 2 \ 

^)b — 

mwsp_c\C^vvii^^ VF, Z) + (^cXyyi{W, Z, W) + s'^C^YVii^^'y^ 

- C24{cz, C^,Cz) - C^C24{c^,Cz, C^) - S^C24(c=^, C^, C=^)] 

- ^mwslsf,_JC^^^,{G^,Z, 7 , P^) + Cj^^AiW, Z, G±) - 7 , G±)] 

- g^ml^^s^_^Co{Z,G^,Z) - g^m^Y^s^0_MW,h,W) - gml,Sf,_^c}_^Co{W,H,W) 


'-W 


+ 9 ‘^-^'i^w^G+G-hCoiG'^, Z, G^) + s^?Ti^Ac+c-/jCo(G^, 7 , G^) 

CvF 

+ 65'^A/i/i/im^s^_^Go(/i, TV, h) + 2g^Xi{}jhrn^(?p_^GQ{H, W, ff) 

+ '^ 9 '^^Hhh'fnw^0-a^P-a[Go{h, W, H) + Gq{H, W, h)] 

3 

+ |-m^s^_„{s2_„[G24(VF, /i, G±) + G24(G±, /i, Py)] 

+ C^-a[<^24(Py, H, G±) + G24(G±, P^) + G24(PP^, K H^) + G24{H^,h, PP/) 

-G24(PF,^^,i^^)-G24(i^^,^, 

+ ^rnw^s^_^[G24{G^, G^, Z) + G24(^, G±, G^)] 


'\G+G-hG24{G^,G^, G±) + XH+H-hG24{H^,A, H^) 

+ 2XGGhG24{G^, G±, GO) + 2A^^;,G24(A, H^,A) 

Q\hhhG 24 {h, G^,h) + 2 XHHhG 24 {H, , H) 

+ ^G+G-hG24{G^, h, G^) + Xu+u-hG24{H^ , H, H^) 


■ 9 


2 2 

■ 5 -5^-a 


- g^cl_^ [QXhhhG24{K H^,h) + 2XHHhG24{H, G^,H) 

+ ^G+G-hG24{G^, H, G^) + Xfj+H-hG24{H^, h, H^) 

- ff'A^+G-h^/3-aC^-a[C'24(G±, h,H^) + G24{H^, K G^) - G24(G±, H^) - G24{H^,H, G±)] 

- 2g^XHhhSp_o.Cp_^[G24{h,G^,H) + G24{H,G^,h) - G24{h,H^,H) - G24{H,H^M 


9 mwSg_a 


3Bo{q^,W,W)+3Boiq^Z, Z) - 4 


9 


„2 


2 ■ ‘G+G-hBo{q^ G±, G±) + ^-XGGhBo{q\G\ G^) + ^XhhhBoiq^ h, h) 

.2 


+ \>^H+H-hBo{q\ + ^XAAhBoiq^A, A) + ^XHHhBo{q^ H, H) 


3 4 

^-:rmwSf,_jBo{pl PT, h) + i?o(pi, PP^, h) + ^[i?o(p?, Z,G^) + Bo{pI Z, G±)] 
+ Sw[Bo{pI, 7,G^) + Bo{pI,j,G^)]Y 


(C54) 
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/ 2 2 n-1t^2,1PI / 2 2 2\ 

[9 O' )b — 


gmwsp-a 


9 2 


CX^^2{Z,W,Z) + c^cX^:^^{w,z,w) + 

- Ci22‘i{cz:C^,Cz) - C^C'i223(c^,C^,C=^) “ 'S^C'l223(c^, , C=^) 


2 

9 

—1 


- 7, + Cl^UW, Z,G^)- 7, G±)] 

+ §rnvys^_„ 2(W", K G±) + J;2(G^, /i, 

wss 

'hVV2 
-iSSV fr^± 


WVS 


+ %mwsp_^cl_^ CWUW, H,G^) + QVy2(W, K H^) - H, H^) 


+ C^^UG^,H, W^) + h, T^) - G^^UH^,H, t^) 

SSV //-yO n± 


+ C]^yy2{Z, G^, G°) + G^yy2{G'^, G^,Z) 


— s 


'^G+G-/iC'i223(G^5 G'°) + ^H+H-hCl223{H^, H^) 

+ 2\GGhCi223{G^, G±, G°) + 2Aa^,,Gi223(71, H^,A) 

^ 6A/i/i/iGi223 {h, G^, /i) + 2\HHhGi223 {H, , H) 

+ -^G+G-/iC'i223(G'^! G^) + ^H+H-hGl223{H^,H, H^) 


— c 


2 

0 —a 


eXhhhGi223{K H^,h)+ 2XHHhGi223{H, G^,H) 


+ AG+G-hGl223(G±, F, G±) + \H+H-hGl223{H^,h, H^) 


- >^H+G-hSf,-aCf3-a[Gl223iG^,h, H^) + Gi223(^^, K G^) - Gi223(G±, H^) " Gi223(^^, G±)] 

- 2\HhhSp_^Cp_AGl223{KG^,H) + Gi223(^,G±,/l) - Gi 223{KH^,H) - G1223{H, , h)]. 

(C55) 

^hWW^P^''P2-!9 )b = 0, (C56) 
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where 


[18(724 + Pi(2(72i + 3(7ii + Cq) + ^ 2 ( 2^22 + C 12 ) + Pi ■ P2(4(723 + 3(7i2 + Cu — 4(7o)] {X, Y, Z) — 3, 

[3(^24 +p\{C 2 l — Co) +P 2 {C 22 — 2(7i2 + Cq) + 2pi ■ P2{C23 — Cu)] {X,Y, Z) — 

CI^UX,Y,Z)^ 

[3(^24 + p\{C 2 i + 4(7ii + 4(7o) + P 2 {C 22 + 2 C 12 ) + 2pi • ^ 2(^23 + 2Ci2 + (7ii + 2(7o)] (X, Y, Z) — —, 
Y, Z) = (IOC 23 + 9 Ci 2 + Cii + 5Co) (X, y, Z), 

Y, Z) = (4Cii - 3 Ci 2 - ^23) (X, y, Z), 

Cl^^^iX, y, y) = (2Cn - 5 Ci 2 - 2Co - C 23 ) (y, y, y), 
y, y) = (6-23 + C 12 + 2Cii + 2Co){x, y, y), 

y, y) = (C 23 - Ci 2 )(y, y, y), 

Ci223(^, = (C 12 + C 23 ){x, y, y). (C57) 

4 . Decay rates for loop induced processes 


The decay rates for the loop induced processes are given by 


T{h 77 ) = 


r(h ^ y7) = 


y/^Cpa, 


em^h 


2567r^ 
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evcX^h 


1 - 
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2 \ 3 


ml 


1287r3 

Sp-aJv + ^ QfN^VfJp - g+g _ s'^)Js 


r(/i gg) = 


y/2GFalmf^ | 

1287r3 1 




The loop functions are dehned as 
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Is = —[1 + 2m^±Co(0,0,m^,mj;^±,mj;^±,mj;^±)], 
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(C58) 
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and 



h '"'W 


Qc^{ml - m|)C'o + 2s^(m| - m|)C'o|, 


(C64) 
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